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We derive and study two different formalisms used for non-equilibrium processes: The coherent-state path 
integral, and an effective, coarse-grained stochastic equation of motion. We first study the coherent-state path 
integral and the corresponding field theory, using the annihilation process A + A —> A as an example. The 
field theory contains counter-intuitive quartic vertices. We show how they can be interpreted in terms of a 
first-passage problem. Reformulating the coherent-state path integral as a stochastic equation of motion, the 
noise generically becomes imaginary. This renders it not only difficult to interpret, but leads to convergence 
problems at finite times. We then show how alternatively an effective coarse-grained stochastic equation of 
motion with real noise can be constructed. The procedure is similar in spirit to the derivation of the mean-field 
approximation for the Ising model, and the ensuing construction of its effective field theory. We finally apply 
our findings to stochastic Manna sandpiles. We show that the coherent-state path integral is inappropriate, or at 
least inconvenient. As an alternative, we derive and solve its mean-field approximation, which we then use to 
construct a coarse-grained stochastic equation of motion with real noise. 


I. INTRODUCTION 

Stochastic processes are ubiquitous in nature: Think of gold 
particles suspended in water, which aggregate upont collision 
111 El, a beautiful realization of the diffusion aggregation (or 
annihilation) process A+A —> A. Think of sand grains rolling 
down a hill, or its cellular automaton representatives, such as 
the Bak-Tang-Wiesenfeld j3] or the Manna sandpile J4J mod¬ 
els. Even simpler, think of a large number of particles diffus¬ 
ing. To understand the physical properties of these systems, 
several routes are open: One may start with a direct numerical 
simulation of say the mentioned gold particles. For techni¬ 
cal reasons this study would be restricted to a relatively small 
number of particles. Thus, in a second step, one strives for a 
more efficient effective description. This could be achieved by 
dividing the system into boxes of size l, counting the number 
of particles inside each box, and trying to derive an effective 
coarse-grained description for the evolution of the number of 
particles inside each box. The question then arises, how do 
we do this? 

Let us step back, and consider an example from equilib¬ 
rium statistical mechanics: In order to understand the phase 
transition between the ferromagnetic and the paramagnetic 
phases in a ferromagnet, or the liquid-gas transition in wa¬ 
ter, one first reduces these phenomena to the simplest possible 
model, in both cases the Ising model. The latter can be stud¬ 
ied numerically, or through analytic techniques. An analytic 
treatment may start from the mean-field approximation, and 
then progress to the construction of a coarse-grained model, 
also termed effective field theory. What one learns from mean- 
field theory enters into the effective field theory as the descrip¬ 
tion inside a box, usually with one or few degrees of freedom 
(fields). This construction has to be supplied with an addi¬ 
tional coupling between boxes, which completes the effective 
field-theory. It can then be analyzed with renormalization- 
group techniques. The latter are expected to give the correct 
universal properties, as e.g. the divergence of the specific heat 
when approaching the critical point, even though the precise 


location of the phase transition temperature itself has been lost 
when constructing the coarse-grained description inside a sin¬ 
gle box. 

Coming back to our discussion of the aggregating gold par¬ 
ticles, the key point is the derivation of an effective field the¬ 
ory. There are two general-purpose methods to do so, both 
with their unique strengths and weaknesses: The coherent- 
state path integral (CSPI), and the coarse-grained stochastic 
equation of motion (CGSEM). In these notes, we will study 
both techniques side by side: 

The coherent-state path integral (CSPI) has proven to be a 
useful tool, both for quantum many-body problems mm, as 
in statistical mechanics CHM- Despite its success, e.g. for 
reaction-diffusion processes, its use led to quite some confu¬ 
sion. Indeed, as we will see below, the CSPI quite naturally 
introduces an imaginary noise, rendering a physical interpre¬ 
tation difficult. The literature on the subject is vast SE3- 
but leaves unanswered key questions the author of these notes 
asked himself. It is his intention to close this gap. 

We start by giving a pedagogical, introduction to the CSPI 
(section III). This is mostly standard, following the work by 
M. Doi GJ[8l, L. Peliti (9), and the beautiful introduction by 
J. Cardy lTT0l . 

For concreteness, we then focus on reaction-diffusion pro¬ 
cesses, as the gold aggregation process discussed in the be¬ 
ginning, and construct a field-theory action. This kind of pro¬ 
cesses leads to the appearance of some surprising, and seem¬ 
ingly counter-intuitive vertices in the field theory, which are a 
consequence of the conservation of probability. We show how 
they can be interpreted in terms of a first-passage problem 
(section|III). 

The field theory can then be reformulated as a stochastic 
equation of motion (section |IV| ). It has an imaginary noise, 
which gave rise to some puzzlement in the literature uni ns- 
usi. As we will show below, contrary to real noise imaginary 
noise leads to a narrowing of the probability distribution. As 
the basis of the CSPI are coherent states, equivalent to Poisson 
distributions, the presence of an imaginary noise tells us that 
over time the probability distribution becomes narrower than a 
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Poissonian distribution. Coding for such narrow distributions 
with Poissonians is only possible via complex states, i.e. inter¬ 
ference. We show, and check numerically, that the stochastic 
evolution of the coherent states allows to sample directly the 
evolution of the discrete probability distribution, starting from 
the initial Poisson distribution 

Pi{n) = ( 1 ) 

n\ 

to the final distribution 

w( ”) = ( e "“'^r) £ ■ (2) 

The average goes over endpoints af of trajectories of the 
stochastic equation of motion, for different realizations of the 
noise £ t . As we will see, this formalism breaks down when (if 
has diffused “too far” away from the positive real axis. This 
is an intrinsic problem of the CSPI, and difficult to repair. 

As mentioned above, there is a second, alternative approach 
(section |Vj. To this aim, one replaces the discrete number of 
particles on a given site by a continuous variable: Either, as 
an ad-hoc procedure, or via coarse graining, introducing the 
particle density in a box of size £. Demanding that the result¬ 
ing continuous random process has the same drift and vari¬ 
ance as the underlying discrete process leads to an effective, 
coarse-grained stochastic equation of motion (CGSEM) with 
drift and real noise. As in the CSPI, its amplitude is propor¬ 
tional to the square root of the drift term, with the difference 
that the latter is real, while the former is imaginary. Contrary 
to the process with imaginary noise in the CSPI, the real pro¬ 
cess converges well for all times, and yields an efficient effec¬ 
tive description. It is at the basis of most effective stochastic 
field theories. However, the stochastic equations of motion 
are rarely derived, even though the procedure given below is 
quite generally applicable. Most often, the stochastic equa¬ 
tions of motion are conjectured on symmetry grounds, more 
obscuring than enlightening their origin. 

We then proceed to a non-trivial example, the stochastic 
Manna sandpile (section m- The rules are simple: If the 
number of grains on a site exceeds one, two grains are redis¬ 
tributed, or toppled , onto randomly chosen neighbours. We 
study this model numerically, and show that coherent states 
are not an appropriate basis for a coarse-grained, stochas¬ 
tic description: On one hand, the probability distribution for 
Manna sandpiles is an exponential, and not a Poissonian, as 
the coherent state. On the other hand, while a Poisson distribu¬ 
tion has one parameter, we observe that Manna sandpiles are 
characterized by two parameters. Thus, while a description 
in terms of a CSPI is always possible (and at least for short 
times exact), it is also plagued by the appearance of complex 
states, and the corresponding convergence problems. In hind¬ 
sight, passing into the complex plane may not be surprising, 
as it can be interpreted as the “trick” of the CSPI to generate a 
second dynamic variable. We then turn to a more efficient de¬ 
scription, and construct an effective stochastic field theory. To 
this aim, we define a variant of the Manna model, the range-r 
Manna model: its toppling rules are modified s.t. grains end up 
not on neighbouring sites, but on sites within a distance of r. 


We show that it converges for r —> oo to a mean field model, 
which we solve analytically. Using this mean-field model as a 
coarse-grained description for an elementary box, we derive a 
stochastic field theory for the Manna model. This field theory 
is known l20lj23ll . The advantage of the present scheme is that 
we do not have to invoke symmetry arguments, and that our 
scheme fixes all parameters, restricting the classes of models 
to be considered to a sub-manifold, which is equivalent to the 
simplest dissipative dynamics of a driven disordered manifold 

M- 

After the conclusion (section [VTll ), the reader will find sev¬ 
eral appendices to which more technical details have been rel¬ 
egated. 


II. THE COHERENT-STATE PATH INTEGRAL (CSPI) 

The coherent-state path integral (CSPI) is a formalism 
which evaluates exactly the evolution of probabilities for a 
stochastic process. To this aim, the different configurations 
of the system are represented as in quantum mechanics by n- 
particle states |n). This allows to write probabilities p(n) as 
states, i.e. superpositions | ijj) := P( n ) l n )- The evolu¬ 

tion operator is then encoded into a Hamiltonian, acting on 
these states. Finally, a path integral is introduced. Its eigen¬ 
states are coherent states, i.e. eigenfunctions of the annihila¬ 
tion operator to be defined below. 

Having constructed an exact representation of the stochastic 
process as a coherent-state path integral, the latter can be stud¬ 
ied with different methods: Either using perturbation theory, 
possibly coupled with renormalization group methods (sec¬ 
tion or by rewriting it as a stochastic equation of motion 
for the states | if) (section |Tvj i. We will study these techniques 
in turn. 


A. Quantization rules 

Consider a single site which can be occupied by n particles 
(bosons), n = 0,1,.... Denote this n-particle state by 

\n) := (at)” |0) , (3) 

where |0) is the normalized vacuum state (0|0) = 1. While 
cfi is the creation operator, its conjugate a is the annihilation 
operator, a |0) = 0. They have canonical commutation rules 

[a, a 1 ] = 1 . (4) 

The scalar product between two states is 

(n\m) = (0| a n (a t ) m |0) = n! S nm . (5) 

This is proven by commuting all a to the right, using Eq. Q. 
Thus | n) is not normalized to 1, but to ( n\ n) = n\. The 
number operator is n := a' a, i.e. 

n | n) = a)a \n) = n\n) . 


( 6 ) 




We note for convenience that 


We thus arrive at 


aa' |n) = (n + 1) | n) , (7) 

a 2 (a t ) 2 |n) = (n + l)(n + 2) |n) , (8) 

(a^) 2 a 2 |n) = n{n — 1) |n). (9) 


B. Master equation and Hamiltonian formalism 

We now want to code a master equation for the occupation 
probability in this formalism. Suppose the probability for hav¬ 
ing n particles at time t is pt(n), with Pt{ n ) = 1- We 

associate with this probability a state 

OO OO 

\A) ■■= J2pt(n) | n) = ^p t (n)(a t ) n |0> . (10) 


n—0 


n—0 


Consider the master-equation for the probability pt(n). 


dtptin) = - (n+ 1 )np t {n + 1) - n(n - 1 )p t (n) 


+P 

+K 


(n + 1 )p t (n + 1) - np t (n) 
{n - l)p t (n - 1) - np t (n) 


( 11 ) 


In the first process two particles meet and annihilate with rate 
v: A + A — A. In the second process, a particle decays with 
rate p: A 0. In the third process a particle “gives birth” to 
two particles with rate k: A —A + A. Note that probability 
is conserved, ]T^ 0 dtPt(n) = 0. 

We now want to derive the “Hamiltonian” associated to 
this master equation, in the form 


d t |'t/’t) = H | ip t ) ■ 


( 12 ) 


To this aim we multiply both sides of Eq. (Ill with (a 1 ) 71 |0), 
and then sum over n. The factors of n are expressed using the 
number operator n = a'a, 

OO 

d t J2Pt(n)^r\0) 

n =0 

oo 

= [ Pt ( n + l )^d 2 a ] - Pt{n)(a^) 2 a 2 (a t ) n |0) 

n—0 

oo 

+ p'^2^Pt{n + l)aa 11 - p t (n)a^a (d t ) n |0) 

n—0 

oo 

^2 \pt{n- l)(a t a- 1) -p t (n)a)d (a t )” |0) . (13) 


+ K 


n—0 


Note that we have taken advantage of relation s f7| > to (|9]> to 
simplify the expression. Next we use definition (TTOJ to rewrite 
this expression in terms of |t/> t ). As an example consider the 
first term on the r.h.s., d^cPptin + l)(a^)" +1 |0) = 

£~ 0 a) d 2 pt{n)(d)) n |0) = a^b 2 |t/> t ). We extended the sum 
to n = 0, which is possible since the first term on the l.h.s. 
does not contribute, due to the preceding operators a 2 . 


d t\A) = 2 “ (a f ) 2 a 2 | i>t) + P 


a — a' a 


I i>t) 


+K 


(a t ) 2 a — a^d \ip t ) ■ 


(14) 


Using Eq. ([T2|, this identifies the Hamiltonian 




'An 2 


-m 


U 2 n 2 


P 


At 


(a^) 2 a — a^b 


05) 

This Hamiltonian is normal-ordered, i.e. all a* stand left of all 
d. It has all the terms expected from quantum mechanics, ex¬ 
cept that for each expected term there is a second term which 
does not change the particle number, and which ensures the 
conservation of probability. Indeed, conservation of probabil¬ 
ity can be written as 


OO 

» = a£ 


•t l^Pt{n) = d t (0|e° | 4h) 

n—0 

t 


= (0| e a H(d' ,a) \ip t ) 

= (0|7f(a t + 1, a) e a \f> t ) . 


(16) 


For the first line we used that the 1/n! in the definition of the 
exponential function cancels the normalization ([5]). For the 
second line we used that 


e Aa /(fi f ) =/(fi f + A)e Aa 
e Aat /(a) = f (a — A)e Aat 


(17) 

(18) 


Noting that an a) inside H , when acting to the left on (0|, gives 
no contribution, we arrive at the constraint of conservation of 
probability for the normal-ordered Hamiltonian 'H 


H(a\a) 


= 0 . 


(19) 


Eq. ( fl9| ) is a necessary condition to ensure that ( fi~6| ) holds; 
it is also sufficient since using (17 1 the state = 

X^ n 7'i( n )(®' + l) ra |0) can be chosen arbitrarily. Looking 
back at Eq. ( [T5| , we see that the second term inside each 
square bracket is such that at a' = 1 the sum of the two terms 
vanishes, thus as stated it ensures the conservation of proba¬ 
bility. 


C. Combinatorics 

Let us remark that the combinatorics used in the above pro¬ 
cesses is the basic combinatorics of choosing k out of n par¬ 
ticles, (™), relevant e.g. for the meeting probability of two 
particles. While this choice is canonic, situations may arise 
where the combinatorics is different. If the stochastic process 
was to contain factors non-polynomial in n, then the Hamil¬ 
tonian would not be as simple, and might e.g. become 
non-analytic in a and a + ; much of the technology developed 
here would no longer work. This holds especially true for the 
stochastic equation of motion to be introduced below, which 
relies on the fact that, via a suitable decoupling, the Hamilto¬ 
nian can be rendered linear in a'. 
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D. Observables 

Now consider an observable 0(n), which depends only on 
the occupation number n. Using the same tricks as in Eq. Unf¬ 
its expectation value can be written as 

OO 

(0)^ t ■= ^ZO{n)p t {n) 

n —0 

= (0| e°0(a t a) | ip t ) 

= (0| 0(a^a + a)e a \ip t ) 

= (0| e> N (a t + l,a)e a \ip t ) 

= (0|O N (l,a)e a |t/» t ) . (20) 

From the second to the third line, we used Eq. {T7]>. In the 
second-to-last line we have introduced the normal-ordered 
version of the operator O, obtained by commuting all a to 
the right and all a' to the left. It is genetically a function of a 
and a', not n = a^a. The last line uses that at acting to the 
left vanishes. 


E. Coherent states 


Coherent states play a key role in the path-integral formal¬ 
ism to be developed below. We define them here, and study 
some of its properties. Coherent states are constructed s.t. 

\<f)) := e^ at |0) =► a\4>) = <j>\<f>) . (21) 


Let us start with 0 real and positive. Then by definition a 
coherent state has a Poisson probability distribution for n-fold 
occupation, 


Pin) = e r . 


( 22 ) 


Note that the definition (21 1 does not contain the factor of e 


thus it is not normalized. This is for convenience reasons; one 
may think of it as a histogram. 

States |0) with a complex 0 are possible too. Since 0 is 
continuous, but the number n an integer, coherent states form 
an over-complete basis, even for <(> > 0. However, not all 
probability distributions can be written as a superposition of 
coherent states with positive weights, i.e. as 


IV>) 



d0/K0)e 0 10) 


(23) 


with p(<f>) > 0. There are several ways out of this dilemma; 
One can use negative (or complex) weights p{(f>), states with 
complex 0, or a combination of both. The formalism to be 
developed below will exploit this freedom. 


By definition, the adjoint state is 


(</>*( = (0 e^*“ . 

(24) 

Eq. (p~7|> implies that the scalar product is 



(25) 


Let us give an interpretation of the adjoint state: Apply (0*| 
to the state \tp t ) defined in Eq. ( |To) >, 

OO 

Pt{4>*) ■= (0* \i ) t) = . (26) 

n—0 


This is nothing but the generating function of the probabilities 

Pt(n), n = 0,1,2,.... 

Now consider the expectation value of a normal-ordered 
observable 0(n) = On (a , a) in a coherent state | <j>). 


(0| e a 0 N (a' 1 ’, d)e^ at |0) 

' rj> ~ (0| e“ lle^ at |0) 

= (0\e^O N ((a^ + l),(a + 4,))e & \0) 
= On(I^) • 


We used Eqs. 0 and ( fl8| , as well as the vanishing of a act¬ 
ing on the vacuum to the right, and a' to the left. To write 
the last line On needs to be normal-ordered. Also note that a 
factor of e :! ' has canceled between numerator and denominator 


of the first line; it is necessary, since our coherent states (21 1 
are not normalized to unity. 

In coherent states, the number of particles is not fixed. We 
show in appendix |B l|that 


e A ™ = e A “ t£ 


= :e 


(e^ — l)(Ta . 


(28) 


The r.h.s. is called normal-ordered and denoted by around 
the operators in question; it is defined by its Taylor expansion 
in a' and a, and then arranging all a' to the left and all a to 
the right, as if they were numbers. E.g. is : (a' a ) 1 : defined to 

be (a°) 2 a 2 . 

Using this relation, or directly the intermediate result 
Eq. ( |B3|l a t o* = 1, and the definition of an observable given 
in Eq. \21) yields 

(e A ")^ = e^- 1 ^ . (29) 

The generating function of connected moments is the loga¬ 
rithm of this function, 

<e Aft >; = (e A - 1)^ . (30) 

This means that the p-th connected moment of the number 
operator n is 


Let us give some explicit examples 

(31) 


= <t> 

(32) 

(« 2 >, 

= ^(1 + 0) 

(33) 

(- 3 ), 

= 0(1 + 30 + 0 2 ) 

(34) 

0 


(35) 


= {h{h- 1)) 0 

(36) 

</> 3 

= (n(n- l)(n-2)) 0 

(37) 
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The last set of relations can also be derived directly, see ap¬ 
pendix [BT] 


F. Many sites 

We now generalize to L sites, denoted i = 1 
with creation and annihilation operators a\ and a, for site i. 
The canonical commutation relations are in generalization of 
Eq. 0 

di,a] 

A state is then encoded as 

OO 

IV 7 ) = E Pt(ni,...,n L )(a\) ni ...(a{) nL | 0 ) . ( 39 ) 

ni,...,riL=0 



Having hopping both from i to j and from j to i with the same 
rate D leads to 

U = (a,j - . ( 45 ) 

Note that by this definition the rate to leave a site in dimension 
d is 2c? x D, and not D. In the continuum limit, and summing 
over all nearest-neighbour sites, this becomes the Hamiltonian 
of diffusion 

■^diffusion = -D / VaJ,Va x . ( 46 ) 

J X 

To avoid overly cumbersome notations, we have set to 1 the 
lattice-cutoff a, which multiplies the lattice diffusion constant 
D by a factor of a 2 ~ d . 

I. Resolution of unity 


We can also construct a coherent state out of single-particle 
coherent states, 

L 

10) ■= <S> l&> ■ ( 4 °) 

i= 1 


G. Coarse-graining 


When constructing effective field theories, one often coarse 
grains, replacing the state variables of several sites by a com¬ 
mon effective variable. For coherent states, this is particularly 
straight-forward: Suppose we have two sites with coherent 
states | (f>i) and | 0 2 ), and we want to know what the probabil¬ 
ity to have n-fold occupation of the combined two sites is. We 
evaluate 


Pcomb(tt) — 'y ' 


n i=0 L 




TTl! 


(n — n\)\ 


— g—(01+02) ('f' 1 


(01 + 020 


( 41 ) 


Thus, combining two coherent states leads to a coherent state 
with the added weights. 


|0l) ® |02) - > |01 + 02) • (42) 


Finally, if we are interested in the probability for n-particle 
occupation of our system of size L given in Eq. ( f40| i, we get a 
state 


|$) 


E& 


( 43 ) 


H. Diffusion 

Consider now the hopping of a particle from site i to site j , 
with diffusion constant (rate) D. The corresponding Hamilto¬ 
nian is (as expected) 

'H = D(a]~a\y i . ( 44 ) 


The path-integral representation we wish to establish is 
based on the coherent states defined in Eq. (21 1 . The key rela¬ 
tion which we are going to prove is the resolution of unity 


1 = — 
27T , 


d0d0*e - ^* |0) (0* 


( 47 ) 


The complex-conjugate pair is 0 = <j) x + i<p y , 0* = 0 X — 
i(j) y \ the integration measure is d0d0* = 2 d(j) x d(j)y. Inserting 
these definitions, the r.h.s. of Eq. <E3 can be rewritten as 

f d0 x d0y —4><h* Acd 


7T 

OO OO 


-EE 

n—0 m —0 


e ^ e^“ |0) (0| e 0 “ 
27r dO roa 


/ o 


7T 


to 


drre~ r r n + m e ie ( n - m ) 

(At)" n m 

X | 0 ) ( 0 | ^ , ( 48 ) 

n\ to! 

where in the last line we set 0 := re . The angular integral is 
vanishing for n/m, resulting in 


OO 


E 

n =0 



d(r 2 ) e 


2\ „-r f 2n 


(«t) 




|0) (0| 


OO 1 

= En( at )"i°) <°i«' 

n\ 

n —0 


( 49 ) 


Applying this expression to the state |to) = ( d') m |0), only 
the term n = m in the sum contributes, and reproduces this 
state. This completes the proof. 

Let us mention another commonly employed trick, namely 
of analytic continuation. This is most prominently employed 
in conformal field theory, see e.g. Ref. If25l . to which we refer 
the reader for the subtleties. The essence is that 0 and 0* do 
not have to be complex conjugates, but that one may think 
of them as two independent variables, which together span 
C = R 2 . A conceptionally convenient choice is 0 real and 0* 
imaginary 1 . 


1 Consider the scalar product (n\ 111 m) ~ f d<pd<p* <)> n (<)>*) m e ^* = 
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J. Evolution operator in the coherent-state formalism, and 
action 

We are now in a position to construct the time evolution in 
the coherent-state formalism. To this aim write the evolution 
operator e 5tn ~ 1 + StH for a small time, and evaluate it 
in the coherent basis, by applying the resolution of unity ( |47| > 
to both sides of e Stn . To avoid problems with not normal- 
ordered terms appearing in (TL) 2 , and higher, we choose St 
infinitesimally small: 

e 8tU{a\a) = ±_J d< p t+std(j) * +st 1_ J d0td ^ 

\ 4 > t+st ) 

x (4>t+8t\ \<j> t ) ■ (50) 


larger to smaller times, we would use the first line of Eq. ( |53| , 
evolving from the final state (0f| to smaller times; the factor 
e -<K+u<l>t+st\ t+St=t{ would then fix (0 f | = (0 t+(5t |. The for¬ 
malism can thus be used both forward and backward in time, 
exchanging the role of 0 and 0*. 

We now consider the forward version, evolving from an ini¬ 
tial state |0i) at t = t\. In the continuous limit, the action from 
time t = t ; to time t = tf becomes 

S[fA) := jf' df tfdttt ~ 0 *] ■ (55) 

(Note that we replaced <f>t +St —> 0| in the Hamiltonian 
'H-Wt+Sf* 0*]’ which is valid in the small-tif limit.) 

The path-integral can then be written as 

|0 f ) = Te^' |0i) 


We need to evaluate the matrix element in question 

(fit+st] e' 5m(at,a) |0 t ) 

~ (0| e^+*‘ a [1 + 6M(tf,a)\ e ^ |0) 

— Q&t + St^t 

x ( 0 | e^‘ at [1 + 5tU{a t + 0 t * +4t) a + 0 t )] e^+*‘ a |0) . 

(51) 


We have used Eqs. 0 and ( [T8] > to commute the exponential 
operators. All operators a are now acting on the vacuum to 
the right, thus do not give a contribution. The same holds true 
for ad acting to the left. Further using the normalization of the 
vacuum state (01 0) = 1, we finally arrive at 


e 5m(at ’ a) |0 t ) = + o{8t 2 ). 

(52) 

Together with Eq. (50 1 , we identify all terms for a time step 
from t to t + St as with 


- St,t+s t 5t = Wt +St - 0?]0t + H(0* +4t , <t>t)6t 

= ti+stlfr - <h+8t] + Htfi+svtMSt. (53) 

The expression S t .t+s L is termed the action for the time step 
from t to t + St. The two possible forms where obtained by 
grouping with either of the factors of e~^ t+st< ^ t + 6t ~^ t ^ t ap¬ 
pearing in Eq. Suppose for the following that we evolve 
from small to large times: Then the second line will be rel¬ 
evant; the unused factor of | t=t . will appear together 

with the initial state 0, as 

e-**# (0| e ^’ a e^ at |0) = , (54) 


J D[ 0 ]D[ 0 *]e- s ^*’^ 



(56) 


Note that (in the simplifying case of one time slice) the state 
10 tf ) corresponds to the state |0 t+4t ) in Eq. (50 1 , thus is part 
of the path integral (i.e. integrated over). On the other hand, 
the state (0j?| in Eq. (50 1 corresponds to (0*. |. When applied 
to |0i), and integrateaover it yields the boundary condition 


0ti = 0i- 

The time-index t at the operators a and «' is introduced for 
book-keeping purposes, to define the time-ordering operator 

T as Te^ dm[a ‘ t ’ 3il := IIt=tf =T putting smaller 

times to the right. (This is the same ordering as in the defini¬ 
tion of the path integral.) 

The final state |0f) is not a coherent state, but the su¬ 
perposition of coherent states |0f). To formalize this bet¬ 
ter, suppose that the initial state is also a superposition of 
coherent states, each with weight p(0i), and normalized s.t. 

POi) : = 2¥ / d( di d0fP(&) = 1. 


W’i) = [ p(</>i) e * |0i) ■ (57) 

J 


Restricting support of p(0;) to 0i > 0 is included as a spe¬ 
cial case, with intuitive physical interpretation. The states 
e _<?ii |0i) are normalized, so that together with the normaliza¬ 
tion of the weight the state |0i) is normalized. Define 


•A(0f|0i) := j P[0]D[0*] e - 5 ^^] 


(58) 


Then 

|0f )=/ |0f) / p(0i) ■ (59) 

J (fif J (f>i 


where we used Eq. {TT}; when integrated over 0*, this identi- B Y construction, |0 f ) is normalized, thus -A(0 f |00 defines the 

fies 0 t as the initial state 0,. Note that when integrating from transition amplitude. 


K. The shift 0* —> 0* + 1 

/d0d0* 0 n (-9 0 ) m e"^* = m!/d0d0*0 n - m e-^*0(n>m). 

For 0 real and 0* purely imaginary the integral f d0*e -<,i ^* yields 5(0), Eq. ( |20[ ) we had considered expectation values of an ob- 

and the former expression vanishes except for n = m. servable O. Suppose we want to measure it at time tf, evolved 
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from \(f>i) at time t,\ until time tf, 

(Otf) = e“* (0|e a C’ N (a t ,a)Te' / ' t i f . (60) 

The factor of e _ ^ ensures that the initial state is normalized. 

Remark now that (0| e a = (1|. Going to the path-integral, 
this can be written as 

(O t{ ) = e“* <1| On^.o )Tett |&> (61) 

fa; =0t 


V[<f>]V[<f>*]0 N ( l,^)e 


4>t-<k P -S[(j>* ,<f>] 


/0 f 


4>t,=4 


The final scalar product yields (1| 0 n( 6 i o) |^f) = 

On( 1, Note that it does not fix (ft = 1, as it would in 

absence of the operator 0^{a', a). 

We now shift all variables <j>* —></>* + 1, to obtain 


- S’ [<!>’’ ,<t>] 


(O tc ) = J V[fi}V[fi*]0 N (lAi)e 
S' [</>*, 4>] := £ dtftdth - W [fit At 


=o 


(62) 

(63) 

(64) 


Note that under this shift 

f 1 d tfidtfa —> 


dt (i tj>* 


= 4> ti -<t>ti+ / dtfitdtfit 


l 


1 )d t (j)t 
u 


(65) 


Apart from the obvious change in the argument of 77, this ac¬ 
counts for the cancelation of the factor of e r?r ~' 9i in Eq. (61 1 . 
John Cardy in his excellent lecture notes ED calls this shift 
the Doi-shift. Its main advantage is that the field tfi* has ex¬ 
pectation zero (at least in the final state). This is particularly 
useful when interpreting the CSPI graphically, as we will see 
in the next section. In addition, the formulae are simpler, and 
more intuitive. Finally, it is advantages when evaluating the 
CSPI via a stochastic equation of motion, see section [TV] To 
distinguish between shifted, and unshifted action, we put a 
prime on the shifted one. In the shifted variables, Eqs. (58 i 
and ([59]) take the form 




4>u—4> f 


V!>i) = 


! 4>f 


e 0f \(j> { ) 


•4'(0f|0*)p(0 i) ■ 


( 66 ) 

(67) 


The initial state is still given by Eq. If one starts from ; 
coherent state then Eq. (67 1 simplifies to 


IV’f) = 


' 4>i 


e ^\fi{} A'(fiflfii) 


( 68 ) 


III. GRAPHICAL INTERPRETATION OF THE 
COHERENT-STATE PATH-INTEGRAL, FIRST-PASSAGE 
PROBABILITIES, AND RENORMALIZATION 


Field theories of the type introduced above are often eval¬ 
uated in perturbation theory, and interpreted graphically. To 


this aim, the part of the action linear in both <f> and <f>* is solved 
explicitly, yielding a single-particle propagator or response 
function. One then starts with n particles, draws their trajec¬ 
tories, and studies how they interact via the terms non-linear 
in <f> and (if. In this process, particles may be destroyed and 
created. In the following, we show how to derive this picture 
from the coherent-state path-integral, based on the shifted for¬ 
mulation in Eqs. ([62|-([68|. 


A. The initial condition 


Start with a general initial state \<f>i, x ) ■= e ^ :c |0). Let 
us create p particles, at positions x \ to x p . In the operator 
picture, this is encoded by 



6 6 





(69) 


Applying the path-integral formalism developed in sections 


IIJ 


and IIK we obtain a field theory with action S ', depend¬ 
ing on the two fields <j> and <jf. In Eq. ( [54] ) we had derived the 
factor for the first time slice, on which we still have to shift 
(f>* -*</>* + 1; this has further to be multiplied by the fac¬ 
tor of e~ Sx4>i,x f rom the normalisation; writing both factors 
explicitly, this results into 




Note the cancelation for the terms proportional to us¬ 
ing Eq. ( [69] ), an initial condition with p particles at positions 
specified above is thus transferred to the path-integral as 




|o) 




(70) 


and 


<k,x = 0 • (71) 

Graphically, we draw a particle emanating from position x at 
time t as a dot at that position in space-time, from which an 
arrow starts, 


1 


X,t 


(72) 


B. The propagator 

Consider diffusion as given by the Hamiltonian in Eq. ( |46| . 
According to Eq. © the action is 

S’o [fi, <F] = [ P x , t dt<f>x,t + . (73) 

J X,t 

This yields the propagator, alias Green, or response function 
(in Fourier space) 

| = = &(t' - t ) e - ^ 4 '-^ 2 . (74) 
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Transforming back to real space, this is 





Q(t' 


t) 


0 4 D(t'-t) 

sj 47t D(t' — t) 

(75) 


It is solution of the partial differential equation 

(d t , - DV 2 ,) Gir! x = S(x - x’)6(t - f) . (76) 

Probability is conserved, i.e. J , G* ,_ x = 1. 


C. The interactions 


To be specific, consider the annihilation process with rate v 
A + A ^4 a . (77) 


The Hamiltonian of this process was derived in Eq. (15 i, 


T-L y [a\a\ = ^ [a 1 ’a 2 - (a t ) 2 a 2 ] . 


The corresponding term in the shifted action is 

w,0] = -E + 

= A E / (C,t + 1) 

x 7t 


(78) 


(79) 


Note that both terms have the same sign. Passing to the con¬ 
tinuum yields 


s'A^A] 



(€,t + 1 ) <&*<!> 


2 

x,t 


(80) 


Since we wish the total number of particles to be 0 X]t —> 
J 4> Xyt , in the discrete version </> xt is the number of particles 
on site x, whereas in the continuum version it is the density of 
particles, resulting in the additional factor of 8 d in the action. 
Alternatively, we could keep o x ^ the number of particles in a 
box of size <5. To avoid these problems, which are not essen¬ 
tial for our discussion, we set <5 —>• 1, except when specified 
otherwise. 


D. Perturbation theory 


Suppose two particles start at time t = 0 at positions X\ 
and X 2 . We want to know the probability pi(tf) to find only 
one of them at time tf. Our formalism gives the following 
perturbative expansion 


Pi (if) = v 





(81) 


The lower two points are fixed at time t = 0, and at positions 
xi and x- 2 - The intermediate times and positions (symbolized 
by black dots) are integrated over. Let us call if the final time. 

Naively, one would expect the probability to be given by a 
path-integral, propagating one particle from X\ at time t = 0, 
and the other particle from x 2 at the same time to a common 
position x at time t, and then propagation of a single particle 
to time tf. Integrating over x, one thus naively expects that 

Pi(t{) = P r d t f G x Zy Gy_ Xl Gy_ X2 

J 0 Jx,y 

= v [ df [ G t y _ Xl G t y _ X2 . (82) 

Jo jy 

This is but the first diagram in Eq. ( [8T] >. The question is, where 
do the remaining terms come from? 


E. Interpretation of the “strange” quartic vertex in terms of a 
first-passage problem 


Let us try to construct the probability of annihilation with¬ 
out making reference to the formalism derived above. This is 
very enlightening, since it will not only re-derive the action 
(791, but also shed light on necessary ultraviolet cutoffs of the 
field theory, and their physical interpretation. This procedure 
is equivalent to renormalization. 

Consider two particles propagating. We draw their posi¬ 
tions Xi(t) for discrete times t = nr, n £ N, represented 
graphically as 



(83) 


Starting at t = 0 at positions Xi(0) = x\ and £ 2 ( 0 ) = £ 2 , 
we want to know the probability pi (r) that the two particles, 
meet, or more precisely are within a distance i5/2 at time t-\ = 
r: 


Pi( T ) 




g:_ x 1 gi,_ x q(\x 


x'\ < 6/2) . 


(84) 


Note that since G(r, x — X\) is not a probability, but a proba¬ 
bility density, we have to say how close they have to come so 
that we consider them to “meet”; the probability that the two 
particles are exactly at the same position is actually zero. With 
this prescription the above expression is a probability density, 
as is G(t. x). To simplify our treatment, we will approximate 
this by (d is the dimension) 



1 2 




























9 


Let us now calculate the probability that the two particles meet 
in the second time step: Particle 1 propagates from x\ at t = 0 
to x' at time t, and then to x at time t = 2r. The second 
particle has intermediate position x". Thus (with the same 
approximation as above) we obtain 



1 2 


( 86 ) 

Now we use that the Green-functions obey the composition 
property 

[ ^X-x'^V-Xl = G 2 J_ Xl . (87) 

J x' 

This allows to rewrite this contribution as 


1 2 

However, this is not the complete result: The particles could 
already have met in the first time step, at time t = r. Sub¬ 
tracting this contribution, we have 




12 12 


Note that one cannot simply subtract the probability to have 
met at t = t. 



12 12 


Let us now calculate the probability to meet for the first time 
at t = 3r, 


Pi(3t) = 


We subtracted the configurations where the particles met at 
times t = t and t = 2r; however this subtracts twice the 
configuration where the particles met at time t = r and at 
time t = 2r, which have to be added at the end. 

Note that once the two particles have met, a single one will 
continue propagating (not drawn here). We can therefore read 
off the action which yields the above perturbation expansion, 



Converting the sum over integer n into an integral, ^ t _ nT —► 
1 fl' d t, and comparing to Eqs. |73| and ( J79j i, we identify the 
rate v as 

v=~ ■ (93) 

r 

It is now clear what the quartic vertex is doing: It converts 
the problem of meeting of the two particles to the problem of 
meeting for the first time, a first passage problem . 

One can try to resum explicitly the perturbation series. As 
we set up the framework, it is well defined for small and 
finite r. Under these cirumstances, resummation is rather te¬ 
dious, and the author of the present notes has decided to elim¬ 
inate the corresponding calculations in order to keep the ma¬ 
terial readable. 

We can, however, deduce the result in the limit of r —> 0, 
and v —> oo: One first realizes that the distance between the 
two particles is again a random walk with a diffusion constant 
2D instead of D. It can thus be described by an action 


5 re i [fi, fi*] 





+ 2DV^ it V4, t ■ 


(94) 


Second, the field &(x, t ) is only defined for x > 0, and zero 
for x = 0: when the two particles meet, a single particle will 
propagate from that point on, and their relative position will 
be zero. This is known as Dirichlet boundary conditions, and 
can be solved with the method of images If26l . In dimension 
d = 1, this leads to 


/-it' — t 

x' ,X 




0(f' 


f*-*') 2 (x+x') 2 

0 8 D(t'-t) — 0 8 D(t'-t) 

a/87t D(t' — t) 

(95) 


Here x is the difference in position at the start, and x 1 the dis¬ 
tance in position at the end. Note the difference to Eq. CD- 
Integrating over x' from zero to infinity, we obtain the proba¬ 
bility that the two particles did not meet up to time t, knowing 
that they started at distance x at time 0, 


Psurvive (*G f) — ^ f 


,V / 8 Dt, 

The probability pi (t) given in Eq. ( |8 T] i then is 

Pi (tf) = 1 ^survive (*£; ff) • 

This can be generalized to higher dimensions. 


(96) 


(97) 


IV. STOCHASTIC EQUATION OF MOTION FOR THE 
COHERENT-STATE PATH INTEGRAL 

A. General formulation 


= J dtdx (/)* t d t fi X! t + D’Ve/)* t V(j) x ,t 
+y X] J + l) ^x,i^x,t • 


(92) 


We established in Eq. ( |66| ) that the transition amplitude be¬ 
tween the coherent states and |</>f) is given by 






(98) 
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Note that we use the shifted action, thus shifted fields <p*, 
since then both <p and (p * have zero expectation values, ren¬ 
dering all following considerations simpler. 

Suppose now that the shifted Hamiltonian, and thus the 
shifted action have only linear and quadratic terms in </>£; a 
term independent of <j>* is absent due to the conservation of 
probability, Eq. ©• 

= #£[&] + ■ (99) 

First consider B[cp t \ = 0, i.e. only a term linear in <p%. Then 
the saddle point obtained by variation w.r.t. <p* gives the ex¬ 
act solution to the path integral, encoded in the equation of 
motion, of (p t 


d t( p t = C[(pt\ , <p ti = fa , (100) 

A'{(pf\(pi) = 6(cpf - <p u ) . (101) 

Quite amazingly, an explicit solution for a non-linear path- 
integral has been given! 

This simple solution is no longer possible if B[<pt] ^ 0. To 
nevertheless use an equation of motion, we introduce a Gaus¬ 
sian random variable, i.e. white noise, £ t , to write 
as an expectation value over the noise, 

(102) 

(ttZt’)t = s(t - t ') . (103) 

Note that if B[<pt\ is negative , then the noise is imaginary. The 
sign of the root is irrelevant, since is statistically invariant 
under Q —► —Q. With the noise, the equation of motion ( |100| > 
changes to 


d t <pt = £[<pt] + VWtlZt , (104) 

(pu = (Pi ■ (105) 

The interpretation is as follows: The transition amplitude 
A'{(pi\(pi) can be sampled by simulating the Langevin equa¬ 
tion ( |104| ), with initial condition ( |105| > and noise ( | 1 ()3[ >, 

A\(pf\(pi) = ^5((p{ - (pt t )^ c ■ (106) 

According to Eq. ( |62| an observable O has then expectation 
at time if 


(0 t{ ) = (0 N (l,(Pt ! )) 




(107) 


This is an intuitive result, with some caveats: First, we re¬ 
mind the replacement of a' —> 1. Second, On (a 1, a) is the 
normal-ordered version of the operator. E.g. isn 2 = ((da) 2 = 
(a() 2 a 2 + (Aa, so that (n 2 f ) = (</> 2 + (pt) £ , see Eq. (33 i. 

In appendix B 3 we give a formal proof of this rela¬ 


tion, based uniquely on the CSPI. The formalism produces 
two terms, a linear term proportional to A^Oi 1, a), and a 
quadratic term proportional to gi- g|~C2(l, a). These two 
terms can then be interpreted as drift and diffusion terms in 


the Ito formalism. This gives an independent derivation of the 
process (|104[ >, and the relation ( |107[ ). 

If \jB[(pt] is real, one may think of equation (104» as de¬ 
scribing what is “ going on ” in the system. This is not the 
case if B[<p t \ is negative, thus y/B[<p t ] purely imaginary: Then 
generically, states sampled by the path integral are “non¬ 
physical” in the sense that they do not correspond to a proba¬ 
bility density, even though the transition amplitude is given by 
Eq. (| 106[>. We will come back to this question in section IV G 


There we will show that in the case of imaginary noise, the 
formalism works for short times, but breaks down for longer 
times. 

In section [V] we will propose a different physically moti¬ 
vated treatment, leading to a coarse-grained effective stochas¬ 
tic equation of motion with a real noise. 


B. Example: Diffusion 


Let us start with simple diffusion, with Hamiltonian 

a] := T~L [d^ + 1 , a] = — D I VbJ,Va x 

J X 


(108) 


Note that the shift has no effect since only VaJ. appears. This 
implies the action, already given in Eq. ( |73j ) 

S' [(p *, (p\ = [ (pl,A<p x ,t + DV<F X)t V<p Xtt . (109) 

J x,t 

Variation w.r.t. <p * leads to the equation of motion 

d t (P x , t = DS7 2 cp Xtt . (110) 

This is a very simple and actually quite remarkable equation: 
While diffusion is a noisy process, leading to fluctuations of 


the number of particles on a given site, Eq. (110 1 is a an ex 


act, noiseless equation. It tells us how the distribution of the 
number of particles on a given site evolves with time. 

As a test, let us check that it keeps the total particle-number 
distribution fixed. Eq. ( |43j ) implies that at, say t. = f;, the total 
particle-number distribution is given by the coherent state 


l*> = 


(111) 


Since for periodic boundary conditions f r d t (p x ,t = 
D f :[ V 2 (p X:t = 0, the state |<E>) does not change over time. 
In particular, this implies particle-number conservation. 


C. Example: Reaction diffusion 


Consider the reaction-diffusion process A + A 
(shifted) action defined by Eqs. (73 i and (|79|>: 


A with 


S'[(P\ </>]=/{ (PlA^t + DW(PlA7(p x ,t 

Jx,t ^ 

+ \ } • ( 112 ) 




















The corresponding equation of motion and noise are 
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dt4>x,t = + DV 2 ^ t + i\Zv(t> X 't£,x,t i (H3) 

(£,x,t£,x',t') = S(t - t')5(x - x') . (114) 

This noise is imaginary. It has puzzled many researchers 
whether this is unavoidable mmmm, or could even 
be beneficial ED- We will come back to this question later. 


D. Dual formulation: Equation of motion for (j>* 


Note that one can define the dual process of Eq. ( |104| >, by 
exchanging in the dynamical action the roles of 0 and </>*: 
Suppose the Hamiltonian can be written in the form 

= £*[#]& + ^b*[#](&) 2 . (ii5) 

The path integral for the generating function at time tf then 
becomes 


M0*):= (e 0 * 0tt ) 

= J T>[0] 


(116) 


Note that this equation is written in terms of the unshifted 
Hamiltonian. Contrary to Eq. ( |6~i~| ) it is normalized, since the 
left-most state is not (0|e a = (1|. Integrating j 0 ^d t 0 t by 
part, and noting that the boundary term changes in the expo¬ 
nential 0*0 tf —► 0£.0i, yields 


p f (0*) = / 1J[(p]V[<f>*]e r ^ i+f t* tdt ** +HKM 




(117) 


This path integral is sampled by the stochastic process 




(118) 

(119) 


It evolves the (dual) state 0^ from tf to t\, backward in time. 


as is suggested by the sign of Eq. (118 i. 

Consider now B* = 0, such that the evolution becomes 
deterministic, —(9 t 0j = £*[ 0 *]- Denote by 'T t ,t { '■ 0 * = 
<J)f —>• 0£, this time evolution, i.e. 




( 120 ) 


Note that 4/ t tf (0) = 0, and \k t4f (1) = 1. 

As a concrete example, consider the branching process, in¬ 
cluding a possible annihilation A —>• 0 


A —^ nA . 


( 121 ) 


Then 


H[a^, a] = X n (a^) n a — a^c 

n 

= f(a*)a - /(l)a f a , 


( 122 ) 


where we defined /( x) := Yh n ^n% n ■ The equation of motion 
( |1 1 8| > then becomes (there is no Doi-Shift) 

- dtti = m) - /(!)</>:, <t>t = r • d23) 


To be explicit, choose A 2 = 1, and all other A* = 0. We have 
to solve (backward in time) the equation dt4>* = 0£ — (0£) 2 . 
It has solution 0£ = 1/[1 — e tf_t (l — 1/0*)]. The function 'T 
then reads 


,t(x) 


x 

x + e tf_t (1 — ic) 


(124) 


Using Eq. ( |117[ ), this yields the generating function, evaluated 
at t = t{, 


p f (0*)=p i (T- tf , ii (0*)). (125) 

This is a classical result, see e.g. 0127]. Suppose one starts 
with a single particle at time t = 0; then pi(0*) = 0*, and the 
above becomes 


Pf (</>*) 


0 * 

0* -f- e* f_ti (1 — 0*) 


(126) 


The probability to have n particles at time tf is given by the 
77 -th series coefficient, namely 


Pn{tf) = e ti_tf (l — e ti_tf ) n \ n> 1. (127) 


This is a rather simple expression. 

We could also try to solve the problem by varying w.r.t. 0*, 
inducing the stochastic equation of motion 


dt<k = 0* + iy/WtZt • (128) 

This equation talks about the evolution of the state |0 t ), who 
will become complex. We will discuss in the next section how 
this can be interpreted. Compared to the latter approach, the 
solution (|126[) is much more elegant, and explicit. 


E. Testing the coherent-state path integral 


Consider the annihilation equation f77j ), 

A + A -24 A 


(129) 


with stochastic equation of motion (1131. For simplicity, we 
concentrate on a single site 2 . 


dt0t = -^0* +*v / W > t& > 
(tef) =6{t-H) . 


(130) 

(131) 


Let us use as initial distribution coherent state |0j), i.e. a Pois¬ 
son distribution with parameter 0;, 


Pi{n) = e 01 -y . 

77 . 


(132) 


There are two ways to study this process. 


2 A complementary study was performed in na for the process A+A ^ 0. 
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P(n) 



FIG. 1: Results for the process (JT29J, using ai = 15, tf — t\ = 0.2, 
and v = 1. Blue diamonds: direct numerical simulation with 5 x 10 5 
samples (partially hidden behind red dots). The blue dashed line is a 
guide for the eye. The statistical error bars are smaller than this line 
width. Red dots: Integrati on of the stochastic equation of motion, 
using SR(Pf S (n)) from Eq. 1 136 1 , for 5 x 10 4 samples, St = 10 -3 . 
The orange dotted lines are %t(Pf (n))± | Q ! ('P f cs (n)) |, defined for all 
real n, which is an estimate of the error. Both simulations ran about 
100s. Within these errors, the agreement is excellent. The dot-dashed 
gray line is a Poisson-distribution with at = 6, which would be the 
result of Eq. ( | 1 30[ l in the absence of noise. Taking into account the 
drift term vat /2 induced by the noise in Eq. ( | 1 29| t (see appendix [c) 
leads to a = 6.43, (black, dotted). The real distribution is centered 
around this value, but is much narrower than a Poisson distribution. 


F. Direct simulation of the reaction process 


Let us start by directly simulating the reaction process 
( |129[ ): First, use the probability distribution ( | 132[ > to obtain an 
integer n (the occupation number at t = tf), and then evolve 
Eq. ( 1 29| ) for a time T = if — ip The latter is best done by 
remarking that if at a given time t there are n particles, the 
probability that they have not decayed up to time t + St (with 
arbitrary St) is 


pT vive (St) = exp (- ” (n 2 . (133) 

Thus one can draw a random number r n £ [0,1], sampling 
the decay probability; solving r n = p™ rvlve (<5t) for St then 
yields 

2 . 

St n ' = 7 yy lnr„ , t n := f Sti , t\ := oo . 

n(n — 1 ) 

v i>n 

(134) 

The t n (with n decreasing) are a series of times when one of 
the n particles decays. Thus the number rif of particles at tf is 
given by 


ttf • tj if 61 tf fj < £ tn f —1 • (135) 

Repeating this procedure, one obtains a histogram for the final 
number of particles, and an associated normalized probability 
distribution pf s (n), where DS stands for “direct simulation”. 



FIG. 2: 5000 samples for the result of the integration of Eq. & 
with v = 1, It — ti = 0.2, (pi = 15. These samples lie approximately 
on a circle of size (j>t = 6.43. For larger times, the samples extend 
further around the circle. 


The result, for (f >j = 15, t{ — t\ = 0.2, and v = 1 is pre¬ 
sented on Fig. |T] (blue diamonds, partially hidden behind the 
red circles). 

Alternatively, one can numerically integrate the master 
equation ( 111 (with p = k = 0) using Pi(n) given by 
Eq. ( 132 1 for n < n max = 3^>, and setting p(n) —> 0 for 
n > u-max- This is the fastest and most precise solution. 


G. Integration of the stochastic equation of motion {Ooj 

Integrating the stochastic equation of motion ( | 130[ > for dif¬ 
ferent noise realizations £ t , and with initial condition & t . = 
(f> i, one obtains a (complex) result for Lf = <j>t r - One then 
measures the final distribution, as an average over all noise 
realizations, 


pf EM (n ) := 


- / 


(136) 


The result is again shown on figure |T] (red circles). One sees 
several things: First, the agreement between the direct numer¬ 
ical simulation of the decay process (blue diamonds, and blue 
dashed line as guide for the eye) and the stochastic equation 
of motion (red dots, and orange dashed lines, with error es¬ 
timate) is quite good. This confirms that Eq. ( 1 36| ) is indeed 
applicable, even though the final states </> f are complex. 

Second, the final distribution is much narrower than a Pois¬ 
son distribution: both the Poissonian obtained for 6 = 6 (gray 
dashed line, result of integrating Eq. ( |130[ ), dropping the noise 
term), plotted on Fig. |T] or the one including a drift term 
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FIG. 3: Result of the integration of Eq. fl30| . with v = 1, t{ — ti = 
0.5, (pi = 15. The black circle has radius <j >f = 3.6614, obtained 
by integrating the drift term dt(j>t = —(fit. + 4>t + t/2. The color 
codes the number of splittings, from yellow over green, cyan, blue, 
magenta to red. (Thus a red point has 2~ 5 times the weight of a 
yellow point.) 


v(fi t /2 in Eq. (130 1 (black dashed line, see appendix [c]for dis¬ 
cussion). Having a distribution narrower than a Poissonian is 
possible only with imaginary noise, which leads to a diffusion 
of the phase of (fit, see figure [2] In contrast, real noise leads 
to a widening of the distribution. We study this in more detail 
below. 

Third, using the stochastic equation of motion has its limits: 
Indeed, already for ti — ti = 0.5, the stochastic equation of 
motion gets appreciable error-bars, even with a large number 
s of samples, and for tf — U = 1 convergence is no longer 
assured. We tried an improved algorithm as follows: Instead 
of starting s “particles” at (fifitfi) = (j>{, and evolving them until 
time tf, whenever one of these particles gets too large a phase 
(which promises to give a large value in we “split” the 

particle in two, each of which carries half of the weight (the 
original weight is w = 1/s) of its “father”. If the phase is 
still too large, we split it again, propagating two particles with 
half the weight each. This procedure is repeated recursively. 
We have not been able to find parameters to improve the pre¬ 
cision at constant execution time. We suspect that when split¬ 
ting points, it becomes more probable that “bad regions” are 
reached, and while the weight of the corresponding points is 
reduced, the probability that they appear is increased. This is 
illustrated on figure [3] We also note that for the toy-model 
studied in appendix C 2 (pure phase) this problem is present. 


This indicates that the convergence problem is severe, and no 
algorithm to overcome it has been found yet. We refer the 
reader to mmm for a more detailed discussion of the 
problems and partially successful attempts at their solution. 


P(n) 



FIG. 4: Result for p GGSEM (n) in Eq. 1 136 1 , after integrating 
Eqs. ( |137[ > (real noise, green circles) and ( | 1 38[ t (imaginary noise, blue 
diamonds) for di = 15, tf — U = 0.025. One sees that real noise 
leads to a broadening of the distribution, while imaginary noise leads 
to a narrowing. 


H. Integrating a stochastic equation of motion with 
multiplicative noise: Real vs. imaginary noise 

Let us integrate the following stochastic equations of mo¬ 
tion: 


d t a t = fi t a t - ^ at , (137) 

d t a t = ifi' t a t + -a t . (138) 

They are constructed such that (|at|) does not change with 
time, see appendix |C| On Fig. [4] one sees that, as expected, 
a real noise leads to a broadening of the distribution (green 
dots), whereas an imaginary noise leads to a narrowing of the 
latter (blue diamonds). For imaginary noise, the phase portrait 
looks similar to Fig. [2] As one can easily observe in a numer¬ 
ical simulation, this leads to problems for larger times, since 
then etf can become imaginary, and both a" and the factor 
of e -af will lead to strong interference between the different 
samples, and to large fluctuations of the such obtained aver¬ 
ages. 

Analytic solutions for both processes are given in appen¬ 
dices |CT| (real noise) and|C 2| (imaginary noise). 


I. Integrating a stochastic equation of motion with 
multiplicative noise: “Canceling” real and imaginary noises 

It is instructive to consider an equation of motion with two 
noises. 


d t a t = fi t a t + ifi' t a t • (139) 

Writing the effective action, and averaging over the noises fi t 
and fi' t yields an exact cancelation of the terms generated in 
the dynamic action: \ J/[a/at] 2 generated from the average 
over fit cancels with — \ f [a/a*] 2 generated from the average 
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n P(n) 




FIG. 5: One trajectory each for process tit, i.e. a direct numerical 
simulation of A + A —¥ A (red, with jumps), and h t , Eq. fl49) (blue- 
grey, continuous, rough). The rate is v ~ 1. We have chosen two 
trajectories which look “similar”. Note that ht is not monotonically 
decreasing. 


over ££. Nevertheless, our equation of motion does not van¬ 
ish. This is possible only since the coherent states are over- 
complete, i.e. we do not need the states | a) with complex a 
to code all possible states. This might allow to define a pro¬ 
jection algorithm, which eliminates the states with complex 
a. 

Let us check this cancelation. To do so write, 

a t = doe** ; (140) 

this implies 

dtOt = ^ + ifr (141) 

for details see appendix[C] Note that there is no drift term (as 
compared to the purely real or purely imaginary cases). The 
probability to find <j> = <f> x + itfiy at time t is given by the 
diffusion propagator, 


e 2t 

P t {(j),(j)*)di(j) x d(t)y = d<f> x d <t>y 

ZTTt 


(142) 


This implies 

P t (a,a*)da x da y = 
1 


2Traa*t 


exp -- 


d(< f> x ,<i> y) 
t'X ; tty) 

ln( —) ln( —)' 
v a 0 ; v a o' 


—- \ P t {4>i f*)da x da y 

o{a x ,a y ) 


2 1 


da x da y . (143) 


The approximation is due to the fact that the larger <p y in 


Eq. (142 1 are not summed over; for pedagogical reasons we 


content ourselves with this approximation. Let us now take 
Eq. (1431, and try to integrate over a. Writing a = (x + iy)a 
a* = (x — iy)a ;, we get 

/ P t (a, a* )e (a ' E+i ° ! ') at daxday 


dxd V c{ x + iy)a.r a t c -l*=i py 2 
2nt 


= e 


(144) 


FIG. 6: Result of a numerical simulation, starting with m = 15 
particles, and evolving for tt — U — 0.025. Blue diamonds: Direct 
numerical simulation of the process A + A —¥ A with rate v = 1. 
Cyan: Distribution of the continuous random walk l |149| i. Red: The 
latter distribution, when rounding m to the nearest integer. Black 
boxes: The size of the boxes in n-direction to obtain the result of 
the direct numerical simulation of the process A + A —> A. Both 
processes have first moment 3.511 ± 0.001. and second connected 
moment 1 ± 0.05; the third connected moments already differ quite 
substantially, 0.75 versus 0.2. 


This result says that real and imaginary noise have canceled, 
thus the superposition of all states gives back the original co¬ 
herent state. This is of course what one expects, knowing that 
the two terms in the effective action cancel. 


V. ALTERNATIVE STOCHASTIC MODELING: AN 
EQUATION OF MOTION WITH REAL NOISE 

A. Stochastic noise as a consequence of the discreteness of the 

states 


In section IV F we had simulated directly the random pro¬ 
cess A + A — v -+ A. Each simulation run gave one possi¬ 
ble realization of the process, in the form of an integer-valued 
monotonically decreasing function n(t). Averaging over these 
runs, one samples the final distribution Pf(n), or, equivalently, 
moments of ri{. Let us now start with a fixed number of par¬ 
ticles instead of a coherent state, n{t\) = rii. We then want to 
ask the question: Is there a continuous random process h(t) 
which has the same statistics as n(t)1 

Let us consider a little more general problem: Be n* the 
number of particles at time t. With rate r + the number of 
particles increases by one, and with rate r_ it decreases by 
one. This implies that after one time step, as long as r±6t are 
small. 


{n t +st - n t ) = (r+ - r-)St , (145) 

((n t+ 6t - n t ) 2 ) = {r++ r-)5t. (146) 

The following continuous random process fit has the same 
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first two moments as nt, 

dh t = (r+ — r_)d t + y/r+ + r_ £ t d t (147) 

(148) 

This procedure can be modified to include higher cumulants 
of n t +st — nt, leading to more complicated noise correlations. 
Results along these lines were obtained in Ref. l29l by con¬ 
sidering cumulants generated in the effective field theory. 


B. Example: The reaction-annihilation process 

In the case of the reaction-annihilation process, the rate 
r+ = 0, and r_ = | fitifit — 1); the latter, in principle, is 
only defined on integer n t , but we will use it for all fi t . Thus 
the best we can do to replace the discrete stochastic process 
with a continuous one is to write 

Mnt - 1) + ^fi t {h t - 1)6 ■ (149) 

Using rii = 15, and v = 1, we have shown two typical trajec¬ 
tories on figure[5] one for the process nt (red, with jumps), and 
one for the process n t (blue-grey, rough). While by construc¬ 
tion both processes have (almost) the same first two moments, 
clearly n t looks different: It is continuous, which n f is not, 
and it can increase in time, which nt can not. One can also 
compare the distribution for if — t[ = 0.5, see figure [6] While 
the distribution of rtf is discrete (blue diamonds), that for n f is 
continuous (cyan). Rounding nt to the nearest integer gives a 
different distribution (red). We have also drawn (black lines) 
the size of the boxes which would produce p(n ) from p(h). 
Clearly, there are differences. On the other hand, it is also 
evident that these differences will diminish when increasing 

Hi. 


C. Diffusion 

We now derive the effective stochastic equation of motion 
for diffusion, i.e. hopping of grains from site i to site i ± 1 
with rate D. This is represented on on Fig. [7] We can not 
directly write an equation of the form ( |147| > for the particle 
number hi on site i, since it does not respect the conservation 
of the number of particles. The latter is realised by introducing 
the current J i+ i t : A positive current J i+ i t = 1 represents 
a particle hopping from site i to i + 1. A negative current 
J i+ 1 t = — 1 corresponds to a particle hopping from site i + 1 
to site i. Each hopping has a rate I)\ the rate for a given 
particle to leave a site is the coordination number 2d times 

D. We thus arrive at the rate equations (with the hat again 
denoting the variables of the continuous process) 

drW.t = dt ( 15 °) 

Ji+ i,t = + \J 1 D(hi : t + hi+ij) Pi+kj 

(151) 



l 


FIG. 7: A coarse-grained lattice with box-size 1 = 4. The yellow 
box contains n = 4 particles. 


The white noise has correlations 

(r7i+i, t ) = 0 , (v i+ i !t V j+ i,t>) = SijS(t-t') . (152) 

To perform the continuum limit, let us introduce the density 
of particles inside a box Bg(x) centered at x and of linear size 
£, as well as the (d-dimensional) current, 

„(*,<):= V J(M):= E hF 

i£Be(x) iGBt(x) 

(153) 

which (in first approximation) is independent of the size of 
the box. (We dropped the hat for convenience of notation.) In 
terms of p, the stochastic equations become 3 (generalized to d 
dimensions) 

d t p(x,t) = -\7J(x,t) , (154) 

J(x,t) = —D\7p(x,t) + \/2Dp(x,t)ff(x,t) , (155) 

(rf(x,t)Tf (x 1 ,t')) = 5 i:i 5 d (x — x')6(t — t') . (156) 

Note that there is no /-dependent factor, neither for the cur¬ 
rent, nor the noise term. Combining the first two equations 
yields 

d t p(x,t) = DX7 2 p(x,t) + 'V[\/2Dp(x,t)ff(x,t)] . (157) 

Let us step back and analyse the above findings; for simplicity 
of notation we again set d = 1. First of all, the diffusion 


3 These equations are standard, undisputed, and appear frequently in the lit- 
erature, see e.g. dj. They are a special case of Eq. (17) of f3H . itself 
equivalent to Eq. (4) of ED. 
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process is constructed such that particles do not interact. A 
given particle will be on a chosen site with probability 1 /L, 
where L is the system size. If N = nL is the total number 
of particles, and n the mean particle number per site, then the 
probability to find n particles on a given site is 



The last relation is valid in the limit of L large. We recuper¬ 
ate our old friend, the normalized coherent state e~^ \<fi), with 
(f> = n. Note that in the CSPI, the analog of Eq. is given 
by Eq. ( |1 10| i, namely 

d t <j>(x,t) = D\7 2 <f>(x,t) . (159) 


Since the CSPI works with coherent states, it does not need 
the noise of Eq. ( | 157[ >. What diffuses is the “weight” (j>(x,t) of 
the coherent state, which is the mean particle number per site, 
termed n above. Thus in the CSPI, both mean and variance 
of the number of grains on a site tends to n. We checked 
with a numerical simulation that Eq. (157 i indeed leads to a 
distribution of grains per site with mean and variance n. 


D. An effective stochastic field theory of the reaction process 


Let us now construct an effective field theory of the 
reaction-diffusion process. Consider a lattice of size L d , with 
particles on it, which can hop from one site to a neighbouring 
one with rate D. To simplify our considerations, let us sup¬ 
pose that we take the limit of v — > oo: if a particle jumps on 
an occupied site, only one of them survives. To construct an 
effective field theory, we introduce boxes of size l. Each of 
these boxes contains n{x, t) particles at time t. 

With rate D, a particle hops. Thus the probability with 
which a particle in a given box will hop is Dn x> tSt\ that it 
will land on an occupied site and thus annihilate is 

ri(x — 1 

5 an n(x,t) ~ D n(x, t)St x ‘ -. (160) 


The second factor is an approximation which neglects the cor¬ 
relations inside the box. If the particle had hopped out of its 
box, then the second factor should involve the density in the 
neighbouring box; writing the density in the same box is an¬ 
other approximation. Last not least, if the box is sufficiently 
large, then one can replace n(x, t) — 1 —> n(x, t). 

To perform the continuum limit, we use the density p(x. t ) 
defined in Eq. ( | 153 i 
then becomes 


The equation of motion of this density 


d t p{x, t) = -Dp(x , t) 2 + \f~Dp{x , f)£(at, t) 

+DX7 2 p(x,t) + V \J2D p(x, t)f}(x, t ) (161) 

(£(x, t)£(x', t')) = S(t — t')5 d (x — x') (162) 

{ri l {x,t)r] 2 {x' ,t’)) = S(t — t')5 d {x — x') (163) 

(£0M)) = (v l (x,t)) = ,t’)) = 0 . (164) 


Note that all factors of St and £ have disappeared, absorbed 
into a non-trivial dimension of £(x, t), and ff(x, t). Note that 
the diffusive noise is usually dropped. The reason is that the 
coarse-grained density p(x, t ) varies smoothly, thus 


V 


y/2D p{x,t)rj{x,t) 


\J2D p(x, t)X/rj{x, t ) (165) 


Integrating the letter over a box of size £ will yield r}{x, t) 
on the boundary, making it less relevant by a factor of l/£. 
It is customarily dropped as subdominant. Thus the effective 
stochastic description of the annihilation-diffusion process is 

d t p(x, t) = —Dp(x, t) 2 + \fDp(x, t)£(x, t)+D'\7 2 p(x, t ) . 

(166) 


E. Comparison between the coherent-state path integral, and 
the effective coarse-grained stochastic equation of motion 


Let us rewrite the equations of motion for the effective field 
theory for the annihilation-diffusion process, first in coherent 
states, and second in the coarse-grained stochastic-equation- 
of-motion formalism, choosing similar conventions. First, the 
stochastic equation of motion for the coherent-state formalism 
reads 

<9t0x,i = + iy/v<l>x,t €x,t + DV 2 (j) x j (167) 

(€x,t) = 0, (£x,t€x',t>) = S(t-t')6 d (x-x') . (168) 

Second, the coarse-grained stochastic equation of motion 
reads (we dropped the noise from diffusion) 

dtPx,t = ~^pt,t + \J~^Px,t£,x,t + D\7 2 p xt (169) 

(£ m > = o , (&,*&;',*> = S(t - t')6 d {x - x') (170) 

Let us summarise our findings, based on what we have done 
so far: 

• Both equations look rather similar. 

• In the formulation with a coherent state <j> x ,t, the noise 
£ X)t is imaginary. This indicates that the distribution 
becomes narrower than a coherent state. 


• In the formulation with the effective density p x t , which 
starts with a sharp distribution, the latter widens by the 
presence of the stochastic noise £ X;t . 

• The noises are both proportional to yjv, the rate, times 
the state variable; they differ by a factor of i\/2. 


• In the coherent-state formulation, perturbation theory 
can be interpreted in terms of particle trajectories, and 
first-meeting probabilities, see section III This inter¬ 
pretation is not possible for the coarse-grained stochas¬ 
tic equation of motion. 
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• The coarse-grained stochastic equation of motion in 
principle has an additional noise term, see the last term 
of Eq. ( fl6l) . This term becomes relevant if we coarse- 
grain with very small boxes, but is irrelevant for large 
boxes: Since it is a total derivative, its integral over a 
volume element is a boundary (surface) term, down by 
a factor of \/l, with £ the box-size. 


F. Other approaches 


As we saw above, the appearance of an imaginary, or 
more generally complex, noise, and its physical interpreta¬ 
tion are puzzling. This is reflected in the literature, see e.g. 
HmEHED- One basic reason is that in the CSPI the vari¬ 
ables are coherent states, i.e. (rather broad) distributions in¬ 
stead of sharp ri-distri hutions. The second reason is that 
they are microscopic, and not coarse-grained variables. We 
then discussed a physically motivated approach based on 
coarse-grained variables. In the field-theoretic context, these 
questions were considered in j29j, and more specifically for 
branching-annihilation in lf32l (see also ed section 9.2). 

An alternative approach is to search a change of variables 
for the coherent states. A beautiful proposal was made in 
Ref. Q3. This approach can be interpreted as rewriting cre¬ 
ation and annihilation operators of the CSPI on each site as 

6* = , a = e~^ p . (171) 


The operator p is the particle-number operator (and not a par¬ 
ticle annihilation operator as a). A derivation of the basic 
equations, and its consequences are discussed in appendix [F] 
As can be seen from the transformation (213 i, if the process 
does not respect particle-number conservation, we obtain ad¬ 
ditional factors of e ±/J ; these terms are non-linear in p ', and 
can not simply be decoupled with the help of a Gaussian noise. 
The situation may be different if particle numbers are con¬ 


served. We come back to this question in section VIF 


Clj 



FIG. 8: Thick lines: The order parameters of the Manna model, as a 
function of n, the average number of grains per site, obtained from 
a numerical simulation of the stochastic Manna model on a grid of 
size 150 x 150 with periodic boundary conditions. We randomly 
update a site for 10 7 iterations, and then update the histogram 500 
times every 10 5 iterations. Plotted are the fraction of sites that are: 
unoccupied (black), singly occupied (blue), double occupied (green), 
triple occupied (yellow), quadruple occupied (orange). The activity 
p = a i(l ~ 1) i s plotted in purple. No data were calculated for 

n < 0 . 5 , where ao = e = 1 — n, Oi = n, and Gq >2 = 0 (inactive 
phase). Note that before the transition, ao = 1 — n and ai = n. The 
transition is at n = n c = 0 . 702 . 

Thin lines: The MF phas e diagram, as given by Eqs. (m) ff. for 
n < |, and by Eqs. (184 1 ff. for n > |. We checked the latter with 
a direct numerical simulation. 


(BTW) sandpile 0. It is defined as follows. 

Manna Model (MM): Randomly throw grains on a lattice. If 
the height at one point is greater or equal to two, then with rate 
1 move two grains from this site to randomly chosen neigh¬ 
bouring sites. Both grains may end up on the same site. 

We start by analysing the phase diagram. We denote by ai 
the fraction of sites with i grains. It satisfies the sum rule 4 


VI. A PHENOMENOLOGICAL DERIVATION OF THE 
STOCHASTIC FIELD THEORY FOR THE MANNA MODEL 


In this section, we apply our considerations to a non-trivial 
example, the stochastic Manna model. We will see that our 
formalism permits a systematic derivation of the effective 
stochastic equations of motion. While the result is known in 
the literature lf20ll23l . it was there derived by symmetry prin¬ 
ciples, which are not always convincing. Furthermore, they 
leave undetermined all coefficients. While many of them can 
be eliminated by rescaling, our derivation will “land” on a par¬ 
ticular line of parameter space, characterised by the absence 


of additional memory terms, see section VID 


A. Basic Definitions 


J2 a i = 1- (172) 

i 

In these variables, the number of grains n per site can be writ¬ 
ten as 



l 


The empty sites are 


e:= a 0 . 

The fraction of active sites is 

a := ^ ai • 

i> 2 


(173) 


(174) 


(175) 


The Manna sandpile was introduced in 1991 by S.S. Manna 

E), as a stochastic version of the Bak-Tang-Wiesenfeld 4 Note that a; has nothing to do with the operator a-i used earlier. 
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We also define the (weighted) activity as 

p:=^2oi(i- 1). (176) 

i> 2 

Note that p satisfies the sum rule 

n- p + e= 1. (177) 

In order to take full advantage of this definition, one may 
change the toppling rules of the Manna model to those of the 

Weighted Manna Model (wMM): If a site contains i > 2 
grains, randomly move these grains to neighbouring sites with 
rate (i — 1). 

On figure [8] (thick lines), we show a numerical simulation 
of the Manna model in a 2-dimensional system of size LxL, 
with L = 150. There is a phase transition at n = n c = 0.702. 
Close to n c , the fraction of doubly occupied sites 02 grows 
linearly with n — n c , and higher occupancy is small. In¬ 
deed, we checked numerically that for n > n c the probability 
Pi to find i grains on a site decays exponentially with i, i.e. 
Pi ~ exp(— a n i), where a n depends on n, see figure [9] This 
is to be contrasted with the initial condition, where we ran¬ 
domly distribute nx L 2 grains on the lattice of size LxL, and 
which yields a Poisson distribution, the coherent state | n), for 
the number of grains on each site, see inset of figure [9] (left). 
This result suggests that coherent states may not be the best 
representation for this system. It further implies that close to 
the transition, p ss a, and we expect that the wMM and the 
original MM have the same critical behaviour. We come back 
to this question below. 


B. MF solution 


In order to make analytical progress, we now study the 
topple-away or Mean Field solution of the stochastic Manna 
sandpile, which we can solve analytically. We define: 

Mean-Field Manna Model (MF-MM): If a site contains two 
or more grains, move these grains to any randomly chosen 
other site of the system. 

The rate equations 5 are, setting for convenience a_ 1 := 0: 


d t CLi = -cij0(i > 2)+a i+2 + 2^^ (aj_i-Oj). (178) 


i>2 


They can be rewritten as 


dtdi = -a,i<d{i > 2) + a i+2 + 2(1 — a 0 — ai)(a,_i - a,) . 

(179) 


We are interested in the steady state dtdi = 0. One can solve 
these equations by introducing a generating function. An al¬ 


ternative solution consists in realising that for i > 2, Eq. (179 1 
admits a steady-state solution of the form 


a,i = a 2 n 


i > 2 . 


(180) 


This reduces the number of independent equations (fa, = 0 


in Eq. (179 1 from infinity to three. Furthermore, there are the 
equations a,; = 1, and idi = n. Thus there are 5 
equations for the 4 variables ao, d\, d 2 , and k. The reason we 
apparently have one redundant equation is due to the fact that 
we already used the normalisation condition \\12\ to go from 
Eq. (fl78| to Eq. ([179 1 . 

These equations have two solutions: For 0 < n < 1, there 
is always the solution for the inactive or absorbing state, 


do = 1 — n , 
di = n , 

di >2 = 0 . 


For n > 1/2, there is a second non-trivial solution: 


do 


di> 0 


1 


1 + 2 n ’ 



4 n 2 — 1 


(181) 

(182) 

(183) 


(184) 

(185) 


(Note that a 2 /d\ has the same geometric progression as 
dj_|_i/dj for i > 2, which we did note suppose in our ansatz.) 
Thus the probability to find i > 0 grains on a site is given by 
the exponential distribution 

p(i) = U _ 1 exp (~ia n ) , a n = In 

(186) 

Using these two solutions, we get the MF phase diagram plot¬ 
ted on figure [8] (thin lines). This has to be compared with 
the simulation of the Manna model on the same figure (thick 
lines). One sees that for n > 2, MF solution and simulation 
are getting almost indistinguishable. We have also checked 
with simulations that the Manna model has a similar expo¬ 
nentially decaying distribution of grains per site, with a decay- 
constant a plotted on the right of figure [9] 

A similar MF analysis can be performed for the weighted 
Manna model, and the Abelian Sandpile Model (ASM); this 
is discussed in appendix [D] 

There is a series of models which interpolates between the 
Manna model, and its MF version: the range-r Manna model, 
where grains are not deposited on a random neighbor, but on 
any site within a distance r. In appendix[E]it is discussed how 
this model converges for large r to the MF Manna model. 


/ 2n + 1 \ 
\2n- 1 ) 


5 After completion of these notes, we found that similar rate equations were 
proposed in the literature for a closely related process, the Conserved 
Threshold Transfer Process, see GE page 229, recursively citing ED, 
ED, and mi We have neither found an application to the Manna model, 
nor the solution (| 1 85^ ff. 


C. The complete effective equations of motion for the Manna 
model 

In this section, we will give the effective equations of mo¬ 
tion for the Manna model. Let us start from the mean-field 
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a 



FIG. 9: Left: (Unnormalized) histogram after manny topplings for n = 2; the probability that a site has i grains decays as e °‘ 585 \ for all 
i > 1. Inset: The initial distribution, a Poissonian . Rig ht: The exponential decay coefficient a as a function of n. The dots are from a numerical 

simulation. The dashed red line is the MF result (214 1 . The green dashed line is a fit corresponding to a — | In ^(n + n c )/{n — n c )^ . Inset: 
blow-up of main plot. 


equations for p(t) and n(t). For simplicity of expressions, 
we use the weighted Manna model. The physics close to the 
transition should not depend on it. Let us start from the hierar¬ 
chy of MF equations ( |D 1 [ ). similar to Eq. ( |179| > for the Manna 
model, and which can be rewritten as 

d t a,i = (1 - i)a.j0(i > 2) + (i + 1 )a i+2 + 2p(a i _ 1 - a,) . 

(187) 

For convenience, let us write explicitly the rate equation for 
the fraction of empty sites e = oq. 


d t e = ci 2 — 2 pe (188) 


The first term, the gain r + = u/> comes from the sites with two 
grains, toppling away, and leaving an empty site. The second 
term, the loss term, is the rate at which one of the toppling 
grains lands on an empty site, r_ = 2 pe. 


We now follow the formalism developed in section V A 
Eqs. ( |145| )-( [T48| . This yields 


d t e = p( 1 - 2e) + y/a ,2 + 2 pe 


t 5 


(189) 


where = S(t — t')/l d , and l is the size of the box which 

we consider. Now remark that close to the transition, a-> « p. 
Inserting this into the above equation, we arrive at 

a t e«p(l-2e) + ^v / TT2i6 , (190) 


Due to Eq. ( | 177[ >. the combination n — p + e = 1, and since n 
is conserved this implies dte = dtp. It is customary to write 
equation ( 190| ) for dtp, instead of d f c. Next we approximate 
\/l + 2e by the value of e at the transition, i.e. e —> e^ F = 2 , 
see the mean-field phase diagram in Fig. [8] We thus arrive at 

dtp- (2n-l)p-2p 2 + JYpit . (191) 


Note that this equation gives back n^ F = and as a conse¬ 
quence of the conservation law n — p + e = 1 also e^ F = 2 , 
used above in the simplification of the noise term. 


Finally, let us suppose we have not a single box of size £, but 
a lattice of boxes, labeled by a d-dimensional label x. Each 
toppling event moves two grains from a site to the neighbour¬ 
ing sites, equivalent to a current 

J(x,t) = —DX7p(x,t) + y/2Dp(x,t)£(x,t) (192) 

with diffusion constant D = 2x^ = 1 The first factor of 

2d d 

2 is due to the fact that two grains topple. The factor of ^ 
is due to the fact that each grain can topple in any of the 2d 
directions, thus the rate D per direction is } >d , resulting into 
D = 1/d. As discussed above, we will drop the noise term as 
subdominant. 

This current changes both the activity p(x, t ), as the num¬ 
ber of grains n(x,t), resulting into d t p(x,t ) = d t n(x,t) = 
—VJ(at, i). It does not couple to the density of empty sites. 
Using the sum-rule \\11\ n — p + e = 1, implies the consis¬ 
tency relation d t p(x, t) = d t n(x,t)+d t e(x,t) for the current; 
this confirms that both p(x,t) and n(x,t) must couple to the 
same current. 

Thus, we finally arrive at the following set of equations: 

d t p(x, t) = ^V 2 p(x,t) + [2n{x,t) - 1 ]p(x,t) 

- 2p(x, t ) 2 + y/2p(x,t) £(x, t ) (193) 

d t n(x,t) = ^\7 2 p(x,t) (194) 

(£(a;, t)£(x', t')) = S d (x — x')S(t — t') . (195) 

This is known as the conserved directed percolation (C-DP) 
class. Instead of writing coupled equations for p(x,t) and 
n(x,t), we can also write coupled equations for e(x,t) and 
p(x,t): 

d t e{x , t) = [l - 2e{x, t)] p(x, t ) + \/2p{x, t) £{x, t) (196) 
d t p(x, t) = ^ V 2 p(x, t ) + d t e(x, t ) (197) 
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The above equations for p and n were obtained in the literature 
112014231 by means of symmetry principles, but never properly 
derived. Evoking symmetry principles also leaves all coeffi¬ 
cients undefined, and does not ensure that Eq. ( |196| > is valid 
on a single site, i.e. is free of spatial derivatives. This locality 
will prove essential in the next section. 


D. Excursion: Mapping to disordered elastic manifolds 


In li24l it had been proposed to use these equations as a ba¬ 
sis for mapping the effective field theory of the Manna model 
derived above onto driven disordered elastic systems. The 
identifications are 


p{x, t) = dtu(x, t) the velocity of the interface (198) 

e(x, t) = F{x, t) the force acting on the interface (199) 


The second equation (197 1 is the time derivative of the equa¬ 
tion of motion of an interface, subject to a random force 

F{x, t). 


d t u(x, t) 


1 

d 


\7 2 u(x, t) + F{x, t ) . 


( 200 ) 


Since p(x, t) is positive, u(x, t ) is for each x monotonously 
increasing. Instead of parameterizing F(x, t) by space x and 
time t, it can be written as a function of space x and interface 
position u(x,t). Setting JF(x,t) —> F(x,u{x,t)), the first 
equation (|196|l becomes 


d t F(x,t) —> d t F(x,u{x,t)) 

= d u F(x,u(u,t))d t u(x,t) 

= l — 2F{x,u{x,t)') d t u(x,t ) 

+ \/2dtu(x, t)£(x, t ) . (201) 


For each x, this equation is equivalent to the Ornstein- 
Uhlenbeck (36) process F(x, u ), defined by 

d u F(x, u) = 1 — 2F(x, u) + V2 £(x, u) , (202) 

(£(x, u)£(x', u')) = S d (x — x')S(u — u') . (203) 


It is a Gaussian Markovian process with mean ( F{x,u )) = 
1 /2, and variance in the steady state of 


< [F(x, u) §] [F(x\ u') - §]) = g 5 V - ^ , )e- 2| “-“' 1 ■ 

(204) 

(205) 


Writing the equation of motion (2001 as 

d t u{x, t) = -\/ 2 u{x, t) + F[x, u{x, t)) 


it can be interpreted as the motion of an interface with position 
u{x, t), subject to a disorder force F{x, u{x, t)). The latter is 
5 -correlated in x direction, and short-ranged correlated in u- 
direction. In other words, this is a disordered elastic manifold 
subject to Random-Field disorder. It can be treated via field 
theory. The latter relies on functional RG (see El for an 


introduction) for the renormalized version of the force-force 
correlator ( |204| i. Functional RG is nowadays well developed, 
and predicts not only a plethora of critical exponents 1581 . but 
also size, velocity and duration distributions (39), as well as 
the shape of avalanches go). 

Also note that Eq. ( |193| l has a quite peculiar symmetry, 
namely the factor of 2 in front of both n(x,t)p(x,t) and 
—p{x,t) 2 . As a consequence, Eq. (1961 does not contain a 
term ~ p 2 {x, t), which would spoil the simple mapping pre¬ 
sented above. The absence of this term can not be induced 
on symmetry arguments only. How this additional term, if 
present, can be treated is discussed in Ref. 


E. Coherent States for the Manna model? 


In the last section, we derived an effective field theory for 
the stochastic Manna model, in the coarse-grained stochastic- 
equation-of-motion formalism (CGSEM). The reader might 
wonder why we did not try to use the CSPI formalism. Well, 
we tried, and here we will share the rather disappointing out¬ 
come: One starts from the discrete Hamiltonian, with rate 1, 
and for simplicity in dimension d = 1 





The first term, proportional to (a| +1 + aj_ 1 ) 2 af checks 
whether there are two or more grains on site i, and then moves 
two grains to randomly chosen neighbours. The last term 
— (aj) 2 a 2 is responsible for the conservation of probability. 
One can show that this is equivalent to the stochastic equation 
of motion 






(207) 


with noise (77 = (€t,j€t>,j') = 8(t - and 

( r lt,j£,t' j') = 0- Note that the first (diffusive) term comes 


from the shift aj —> aj + 1 , to be done before decoupling the 
path-integral with noise terms. 

This equation is quite ugly: Not only does it have two mul¬ 
tiplicative noises, one of them is imaginary, inducing the con¬ 


vergence problems mentioned in section IV G 


As we do not know how to proceed from here, let us try to 
address this question on a more abstract level: Could coher¬ 
ent states indeed help us? Let us also state our prerogative, 
namely that one wants to derive an effective field theory for a 
local field. First of all, one has to realise that what the coherent 
state does is to re-write the number of grains on a given site as 
a superposition of Poisson distributions (aka coherent states) 
on this site. This complicates matters; after all, we want an 
effective variable which gives us some intuition. 

Let us take a step back and look at the derivation of the 
field theory of the Ising model. One realises that there one 
needs coarse graining. What we have done in the preceding 
section was to construct coarse-grained effective variables for 
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a block of t x l sites. As for Ising, this necessitates to make 
approximations inside a block. We did this, supposing that 
we are close to the transition, and that we can use the MF 
approximation for the block. 

But maybe this block can be described by a coherent state? 
From Fig. [9] we see that the distribution inside a block is, at 
least approximately, an exponential, represented by (we nor¬ 
malised) 

i° 0 ) /?>= a o io)+ (1 ~ 1 °! ) ^r /?) at i°) • ™ 

This implies n = y-p , and p = in the example of 

Fig. [9j a 0 « 0.195, /3 « e" 0 - 585 , n = 2, and p « 1.115. 
This has to be contrasted to a coherent state, with on average 
h particles, (also normalised) 

|n) = |0) . (209) 


This state has ao = e~ n , possesses n particles, and activ¬ 
ity p = n + e~ n — 1. There are two problems: First of all, 
the tail is different; this should not be an issue close to the 
transition, where triple and higher occupancy are unimpor¬ 
tant. The second and more fundamental problem is that while 
the state ( 208[ > has two independent parameters, the coherent 
state ( |209 > has only one! Still, the coherent-state functional 
integral will correctly propagate a state. More precisely, it will 
calculate a transition probability from an initial state to a final 
state, after decomposition into a coherent-state representation. 
It will do so by passing through complex intermediate states. 
The imaginary part of these complex variables provides the 
second, “missing” variable. Thus thinking of a single (real) 
coherent state per site is not appropriate. 

One could also try to work with coherent states for the num¬ 
ber of empty, once, twice, or triple, ... occupied sites inside a 
box. This would yield a state of the form 


l a o 


ai, a 2 ,...) := e * 2 (“oa++a 1 a++a 2 at+...) ^ 


( 210 ) 


where the a\ create an i times occupied site. If = 1, 

the expectation of the number of sites will still be £ 2 , but it 
will be a fluctuating variable, with variance L While this is 
probably acceptable, the author does not see what would be 
gained in terms of simplicity of derivation, knowing that all 
the approximations necessary for the MF treatment of section 
IVIBI would have to be made as well. 

Coherent states as a basis for a stochastic description of 
the Manna model have also been proposed by Pastor-Satorras 
j20l , based on the field theory of Wijland-Oerding-Hilhorst 
ED. The idea there was to introduce two species A and B, 
where B particles are associated to the activity, and can dif¬ 
fuse with rate D, whereas A particles are stationary. The par¬ 
ticle number in the Manna model is associated to the total 
number of A and B particles. The rate equations proposed to 

mimic the Manna model are B —4 A, and A+B —A B + B. 
This yields a Hamiltonian 


H[a\ a, b\b] 


+ h -fit) b 
(tf) 2 ba-tfa r ba] . (211) 


+&2 


This theory was then analyzed in terms of shifted fields, see 
Eq. (1.15) of ED, which somehow obscures the analysis. 
What was not realized is that the proposed stochastic inter¬ 
pretation necessitates an imaginary noise. This can be seen 
from the fact that the terms quadratic in a' and b' are of the 
form 



+ linear terms in a', b' . (212) 

The matrix has both a positive eigenvalue corresponding to a 
real noise as well as a negative eigenvalue corresponding to an 
imaginary noise, as is the case for Eq. ( |207[ ). From our above 
considerations this is not surprising. It shows once more that 
the CSPI formalism does not yield stochastic equations of mo¬ 
tion with a purely real noise. 

Finally, let us mention another peculiarity of the Manna- 
Hamiltonian in Eq. ( |206[ ): The rate for a site with n > 2 
grains to topple is proportional to n(n — 1); this is dictated by 
the demand to have a Hamiltonian as simple as possible. The 
original Manna model has a constant rate, while the weighted 
Manna model we defined above has a rate proportional to 
n — 1. While this should be no problem at the transition, it 
imposes a specific rate not present in the original formulation, 
a rate found in pairing (or meeting) probabilities of n particles 
in a box. The latter is indeed the framework in which the com¬ 
binatorics of the coherent-state path integral is appropriate. 


F. A canonical transformation for the Manna model? 


Some efforts were spend to find exact representations of a 
stochastic process, without using the CSPI formalism. The 
hope was that such a reformulation could be interpreted as a 
stochastic process with real noise. In section |VF| we discuss, 
and in appendix[F]we rederive the proposal of Ref. IT71 . This 
approach can be interpreted as rewriting creation and annihi¬ 
lation operators of the CSPI on each site as 



Applying this transformation to the Manna Hamiltonian (2061 
yields 


njP 

'‘'Manna L 


[p\p\ = \f +e^-rt) 2 p l (Pi - 1) ■ 


(214) 

Taking the continuum limit, and dropping higher-order terms 
in the lattice cutoff, one arrives at 


n 


p,c ont 
Manna 


[p\p\ 


f x (Vp\x)) 2 +V 2 p\x) 


p{x)\p(x)-l] . 
(215) 


The action 5^ anna = f xt p*(x, t)d t p{x, t ) -7i[p*,p] allows 
for an interpretation as a stochastic equation of motion with a 
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real noise 77 * ( x, t), 

d t p{x, t) ~ V 2 [p(x,t)(p(x,t) - 1)] 


a/2V 


fj{x, t) p{x,t)(p{x,t) - 1 ) 


(rf(x, t)rf {x',t')} = S 1;i 6 d (x — x')5(t — t') 


(216) 

(217) 


This equation is very similar to the linear diffusion equation 


(157 1 , except that on the r.h.s. the particle number p has been 


replaced by p(p — 1). 

The question remains whether passing from Eq. \2\A\ to 
(215 1 is justified. In the path integral, p is the number of par¬ 
ticles. As a discrete number, p strongly fluctuates between 
nearest neighbors, and one expects jr to do the same. Thus 
the approximation from Eq. ( |214[ ) to ( |2 1 5| > is probably not 
justified. One should first construct coarse-grained variables, 
which would probably lead to a stochastic equation of motion 


different from Eq. (2161. We leave exploitation of these ideas 
for future research. 


VII. CONCLUSION 

In this article, we started with the coherent-state path inte¬ 
gral (CSPI) for stochastic systems, which we then reformu¬ 
lated as a stochastic equation of motion. We showed how the 
evolution of the probability distribution can be followed, de¬ 
spite the appearance of imaginary noise. Limitations of this 
formalism were discussed, especially its (at least practical) 
breakdown at finite times. We also showed how some of the 
appearing vertices can be interpreted as transforming a simple 
diffusion probability into a first-meeting probability. 

We then constructed a complementary formalism, based 
on an effective coarse-grained stochastic equation of motion 
(CGSEM) for a continuous variable. Demanding that drift 
and variance for the underlying discrete system are correctly 
reproduced by the CGSEM fixes the latter continuous process 
uniquely. 

We should stress again that while both the CSPI and the 
CGSEM formalism share some common features, they should 
not be confounded: It is tempting to derive stochastic equa¬ 
tions of motion in the CSPI formalism, and then to interpret 
the coherent state \(j>), i.e. a Poisson-distribution with expec¬ 
tation 4>, by the state <f> itself, equivalent to a ^-distribution at 
(/>; as is used in the CGSEM formalism 11201 . We remarked on 
the example of the reaction process A + A —► A, that starting 
from a Poisson distribution, the probability distribution be¬ 
comes narrower than a Poissonian, which in the CSPI is only 
possible with complex coherent states, thus imaginary noise 
in the equation of motion. On the other hand, the probability 
distribution becomes broader than a (5-distribution, necessi¬ 
tating a real noise in the CGSEM. Both noises have, up to a 
factor of \pl the same strength. 

We concluded our considerations by analysing the stochas¬ 
tic Manna model, and gave a straightforward derivation of 
its effective stochastic field theory. Our procedure, based 
on coarse-graining, fixes all amplitudes, including the noise 
strength. Compared to earlier derivations, this derivation is 


simple and transparent; that it fixes all constants is an addi¬ 
tional advantage. It also ensures the simplest mapping on dis¬ 
ordered elastic manifolds. 

We hope that our work helps to clarify the origin and inter¬ 
pretation of stochastic field theories, and that the techniques 
presented here are more broadly useful. 
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Appendix A: Fermionic coherent-state path integral 

The coherent state path integral can also be applied to 
fermionic degrees of freedom, i.e. to states which can only 
be unoccupied, or simply occupied. This can be done via the 
coherent state path integral with Grassmann numbers. Here 
we only write down the basic equations. 

We introduce fermion creation and annihilation operators 
with 


{c, c f } = 1, {c,c} = {c t ,c t } = 0 . 

(Al) 

The states are 


|l)=c t |0), c 0) = 0 . 

(A2) 

This implies 


(l|l) = (0|cct|0) = (0|{c,ct}|0)= (0| 0) = 1 

• (A3) 

Coherent states 


\ijj) = e^ ct 0) = (1 + ipj) 0) 

(A4) 

(t/>| = (0 e c ^* = (0| (1 + cip*) 

(A5) 

(ip*\ip) = 1 + ip* ip ■ 

(A6) 

Resolution of unity 


J dip*dip | ip) {ip*\e~ r ^ = 11. 

(A7) 


This relation can be checked by applying it to 10) and to ip' 11) . 


Appendix B: Proofs of some relations used in the main text 
1. Proof of Eq. ((28} 


Eq. ( [28] ) reads 

gA n — ^AcPd . ^(e A — l)(Pd. 


(Bl) 
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The proof of this equation consists of two steps: Applying the 
l.h.s. of Eq. (|B1[) to a coherent state yields 


e Aata e 0at |q^ = ^ 
n—0 


xat a 


| 0 > 


= g e A„^ t L |0) 

' n\ 

n=0 


= e e ^ at |0) . 


Thus 


((/)* | e Aat “ |0) = (0| e ?i * a e e ^ at |0) = e e>< ^ . 


(B2) 

(B3) 


On the other hand, 

(0*| :e (e "" 1)ata : |0> 

OO 

= ]T(0|e^ a (e A - 

n=0 

oo 

= ]T(0|e^ a (e A - 

n=0 


i)’ 


(a f ) n a n 


e^ a |0) 


1)" |Q) =e e ^*^ . 


(B4) 


where the expectation is as in Eq. ( |B6[ >. To simplify the 
calculations, we show that if 0[a\a] is normal-ordered, 
0[a\a} = C9n[«' , a], then 0[a\a ] can be replaced by 
0[1, a]: indeed, Eq. ( |B8[ > is proportional to 

(0| e“ (0(a\ a)H[a\ a] — 'H[a\d\0(d\a)) ... 

= (0| e a (0(1, a)U[a\a} - U[a\a]0( 1, a))... (B9) 

since (0| e a O(a t ,a) = <0| e a O(l, a), and (0| e h U[a\a) = 0. 

Thus 


_d_ 

dff 


(Ot t ) 


0(l,a),H[a\a] 


(BIO) 


Next, the commutator can be calculated by remarking that 


0(1, a), / H[a\ d] 


= w(o(l,a),H[a t ,a]) . 


(Bll) 


The operator W denotes all possible Wick contractions be¬ 
tween 0(1, a) and 'H[cv : a]. To proceed we suppose that the 
Hamiltonian has only a linear and quadratic term in a', i.e. 


This proves relation (B1». 



(B12) 


2. A different derivation of Eq. p5) 


with B xy [a] = B yx [a]. Then the commutator will be 


The set of relations ( [35] ) ff. can also be derived from 

V = 0 P (0| e a 10) = (0| e a a p e^ at |0) e ~+ 

= (0| (a* + l) p e a a p e^ at |0) e“* 

= <0| e a (a t ) p a p e^ >a+ |0) e~* 

= (0| e a (h — p + 1)... (h — l)ne^ nt |0) e“^ . (B5) 

3. Formal derivation of the evolution of expectation values in 
the coherent-state path-integral 

Consider the expectation value ( [60] > at time if 

(0^) = e _0i (0|e a O(a t ,a)Te / ‘i f dm[a J’ a * ] |0.) . (B6) 


0(l 1 a) 1 'H{a\a] 


S w[a ’ a| i 0(1 ' a) 




6 S r „ t 6 6 

5a\ Sa'y $ a x Sa y 


0(1,a) 


(B13) 


As a consequence, the expectation (BIO i evaluates to 


it, <£> " > “ \l ( £ ” |al + i®'-" |al ) 4a I ° (1,a) 

(B14) 


1 

+ 2 


We can give an interpretation in terms of a stochastic process 
ft defined by 


We are interested in its temporal evolution. At a slightly 
smaller time, the observable O had expectation 

(Oust) (B7) 

= e - * (0\e a e 5m ^^0(a\a)Te ft T Stdtn ^’^ |0;) . 


d t f x ,t = C x [f] + f B x Jf] + £ Xtt (B15) 

(£x,t£x',t') = 5(t~t')BxA4>\ (B16) 

Of f = 0(1,0 tf ) . (B17) 


Note that we have been able to add the factor Q St ' H W,a\ = 
e <5fH[a] f ,a tf ] s j nce y[af,a\, when applied to the left to (0| e a 
vanishes, see section |HBl Thus the time derivative of the ex¬ 
pectation value of an operator is given by its commutator with 


d 

dff 


(O u ) 


0(a\a), %[ 



(B8) 


Indeed, applying the Ito formalism (see e.g. 11421 ) to the ex¬ 
pectation (O tf )j yields Eq. (B14 1 . Comparing Eq. <B15 i to 
Eq. ( |104[ ), we note that the drilt term also contains a term lin¬ 
ear in B. The reason is that to arrive at Eq. ( |104[ ) the shifted 
Hamiltonian ha d been used. Indeed, this is accounted for by 
shifting in Eq. (B12) ed —>■ ed + 1, C x —> C x + J &x,y, 
while B x , y remains unchanged. Thus Eqs. (104» and (|B 15 1 
are equivalent. 
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Appendix C: Stochastic equations of motion with purely 
multiplicative noise 

The terms in the action which cause us trouble are of the 
form 

8S[a\a}= l -j{a*Ya 2 t (Cl) 

Decoupling with a noise £ t , with (£ t ) = 0, (itCt) = 8(t — t') 
we have 

e -tfS[a*,a] __ (£2) 

Thus the equation of motion will have a term of the form 


d t a t = iatit + •■• (C3) 

This is a multiplicative noise. In the following, we will solve 
this equation, dropping all further terms, both for a real noise, 
and for an imaginary noise. To reserve the notation p t for the 
phase, we noted this variable a*; it reminds the coherent-state 
variable a. 


1. Integrating a stochastic equation of motion with purely 
multiplicative noise: Real case 

To start, consider the easier case of a real noise 


The probability is obtained by inverse-Laplace transforming, 

P(ln(a/a 0 ), t) = j ^ ^ e iA[ln(a/a 0 )-ln(a t /a 0 )]^ 


dA 


i\ ln(a/ao) — A(A+i)t/2 


27r 

[ln(a/a n ) + t/2] 2 


\/2nt 


(CIO) 


This leads to the probability as a function of a, 

[ln(a/a Q ) + t/2] 2 


Integrating f a n P(a, t) we find 

(a") = <e Il(n - 1)t / 2 . 
This yields, as it should. 


'j2 , Kta 


dt <a?) - 


i(n — 1) 


’>• 


(Cll) 


(C12) 


(C13) 


2. Integrating a stochastic equation of motion with purely 
multiplicative noise: Imaginary case 


d t a t = itAt ■ (C4) 

We can make the ansatz 

at = «oe 0t ■ (C5) 

This leads to 

dtpt = it - i • (C6) 

Note the drift term; using Ito calculus with d B t = it.dt, the 
latter equation reads 

dpt = d B t - idf . (Cl) 

We infer, as claimed. 


For an imaginary noise as in Eq. (C3 1 , we start from 


d t a t = i£ta t . 


We make the ansatz 


at — aoe 




(C14) 


(C15) 


This leads to 


dtpt = it • (C16) 

Note the drift term, analogous to the one in Eq. ( |C7| >, but with 
opposite sign. The probability to find p is then given by 


da t = aode 


= aoe 0t 

dpt + ^ dpt + ■■■ 

= a o e 0t 

d B t - ^d t + id B 

= a t dB t 



(C8) 


Therefore, for arbitrary A, the generating function is 

= ( e A ^> = ^ e A 7o dl -tfr-!)^ =e A(A-i)t/2 


(C9) 


P ^ t) = e m- (C17) 

The generating function for at, and its complex conjugate a%, 
or more precisely for their logarithms, reads 


Z(A,A*) = (af(an A } = (a 0 e t / 2 ) 


A+A* 


3 —(A+A*) 2 t/2 


(C18) 

In section IV G we had seen that the probability for n-particle 
occupation is still given by the normalized coherent state. 
Even though here we are only dealing with a toy version of 
the real equation of motion, namely Eq. (C14 1 , we can still 
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FIG. 10: The pro babil ity JC24l (red dashed), compared to a direct 
simulation of Eq. i C14 1 with 10 samples, St = 10 -4 , t = 2, ao = 2 
(blue); the thin vertical red lines delimit the domain of the analytic 
solution. Some simulation points lie outside, due to a finite St. 


For small times, it is dominated by the term n = 0, 


_ acos(b) 2 

P?(b, t) « P b ’ app (6, t) := , . (C23) 

r r V 1 - b 2 ) 


Note that 


e - t /2 f a e -*/2 \ 

P" (Or, t)=~ - Pr - -, A . (C24) 

ci 0 \ a 0 J 


We checked this result numerically, see figure [TO] 


Appendix D: MF solutions for the weighted Manna model, and 
the Abelian sandpile model 

1. MF solution for the weighted Manna Model 

The MF rate equations for the weighted Manna Model are 


study the same observable. This will shed light on the conver¬ 


gence issues noted in section IV G 


M*) = (e-°*-V 

n\ 


= £(-0 


k =0 
oo 





— (n+fc) 2 t/2 


(C19) 


k -0 


For large times t, this converges towards po(t) —> 1 — ao, 
Pi(t) —► ao, Pn>2 -A 0, which is no longer a probability for 
ao > 1. Another option is to write the above expectation as 
an integral over P(<j), t), 


Pn(t) 



e~ a ° e 




(a 0 e ict,+t / 2 ) n 

n\ 


(C20) 


When the time t becomes large, the integral starts to oscillate. 
As an example, for ao = 2 it can no longer be done for tft 5. 

Using 4Cl7> , we can also study the probability distribution 
for the real part 


a r 


a T a* 
2 


Defining 6 := e 4 / 2 a r /ao, this is 


(C21) 


dtdi = -(* - l)a l 0(i > 2) + (* + l)a i+2 

+ 2 [£ay(i- 1 ) (ai-i-cbi). (Dl) 


3> 2 


To give an analytic solution of this set of equations is left as a 
challenge to the reader. On the left of figure [IT] we show the 
result of a numerical integration of the rate equations, starting 
from a Poisson distribution with on average n grains per site. 
As can be seen from the right of figure [TT] the tail of the en¬ 
suing distribution remarkably no longer is an exponential, but 
close to a Poissonian. This can be seen from Eq. <(DT}: Setting 


at —> 


a 2 {2pY 2 


(D2) 


the leading term in i is canceled. Indeed, the numerical solu¬ 
tion confirms that this is approximately true; more precisely, 
(i+i)a- +1 1 f° r * d n practice, one can get the first 

20...50 coefficients, depending on n.) 


2. MF solution for the 2D-ASM model 

As a further excursion, consider the 2-dimensional Abelian 
Sandpile Model (ASM), which moves four grains to the four 
nearest neighbours, if the height at a given sites reaches or 
exceeds four. Its MF rate equations are 


P T b (b,t) = 


1 


£ 


y/2nt(l — b 2 ) 

[acos ( b ) + 2 7r n] ^ 


+ e 


[acos(b) — 27 ttj,]^ 


03 ( a ^ fel ,e- t / 2 ) +03 e -‘/ 2 

y / 27T<(l — b 2 ) 


d t a i = -a i Q(i>4:) + a i+4 + 4^2a j (a^i-aj). (D3) 

j> 4 


Using Y2i a i = 1> one can eliminate the infinite sum in the 
square bracket. Let us again suppose geometric progression, 
this time for i > 4, 


CLi = a^K 


i —4 


(C22) 


i > 4. 


(D4) 
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FIG. 11: Left: MF phase diagram of the Weighted Manna Model (wMM) (thick lines); compared to the MF phase diagram of the standard 
Manna Model (thin lines). Right: The distribution of a, for the wMM as a function of i (red points), for n = 0.8. The blue dashed line and 
green points are a Poisson distribution with 2.286e -0 ' 989l /i!. 


Solving as for the Manna model the 5 first equations, as well 
as the constraints \\12\ and ( |173| ), we obtain a non-trivial so¬ 
lution, with non-negative coefficients, for n > |, 


ao 


a i 


a 2 


at 


1 

An — 2 
2(n — 1) 

(1 - 2 n ) 2 
12(n — 2)n + 13 
2(2n- l) 3 

4(n — l)[4(n — 2)n + 5] [271-3 
(1 - 2 n) 4 [2n- 1 


(D5) 

(D6) 

(D7) 


i> 3 . 
(D8) 



This is depicted on figure [l2] For " < §. any set of positive 
a,, s.t. cii = 0 for i > 3, and which satisfy the constraints 
E» °i = 1 and Ei ia,i = 7i is possible. Which solution is 
picked is given by the initial conditions. Integrating the rate 
equations ( |D3[ ) starting from a Poissonian distribution with 
expectation n leads to the solid lines drawn on figure [12] This 
reproduces the solution ( |D5| >-( |D8] > for n > n c ~ 1.542. We 
have checked these analytical predictions by a direct numeri¬ 
cal solution, see dots on figure [l2| 

Note that the same analysis can be done for an ASM on 
a lattice of coordination number 71 ; the most interesting such 
case is n = 3, e.g. the 2-dimensional honey-comb lattice. The 
MF equations of motion then read 


d t aj = -ai®{i > n) + ffli +Tt + n[y~[aj (aj_i-a»). (D9) 

j>n 


Appendix E: The topple-away Manna model, its convergence to 
MF, and limit of large dimensions 

The question arises, how one can go from the standard 
Manna model with nearest-neighbour topplings to its mean- 


FIG. 12: Main plot: The phase diagram for the MF version of the 
2d-ASM. The branch starting at n = 1.5, and ai ^<3 = | is the MF 
solution m ff. The branch starting at n = 0 was obtained by a 
direct numerical integration of the rate eqs. setting to zero all 
cti,i>ioo, and waiting until a steady state is reached. The dots are 
the result of a direct MC simulation of the MF-ASM model, starting 
from a Poisson distribution; we randomly distributed n times size 2 
grains on a lattice of size 150. Note that Eqs. ( |D3| l allow to predict 
the state reached by Monte Carlo. Also note that the branch with 
at : i >4 = 0, which exists up to n < 3, remains attractive beyond 
n = | (see inset, solid lines). The dashed lines are the solution 
(D5}ff. 


field variant. There is indeed a series of models, which does 
this interpolation. To this aim, define 

The topple-away Manna Model of range r (r-MM): If on a 

site (i.j) there are two grains, do twice: Randomly pick a site 
( i',j') with |i — i '| < r, | j — j'\ < r, and move a grain from 
site (i,j) to site (We exclude the origin; two grains 

may end up on the same site.) 

We can show via a numerical simulation that this variant 
converges to the MF-Manna model for r -T oc; in practice 
for r = 15, one is already very close. On figure [13] we show 
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FIG. 13: log-log plot of n c — 0.5 as a function of the range r. The 

shown line is a power-law fit n c - 0.5 = 0.168069/r. “ FIG - 14: The P hase diagram for r = 6 . The transition is at n c = 

0.517. Figure to be compared to Fig. [ 8 ] The thin underlying lines are 
the MF solution, the thin dotted lines the unstable solution. 


how n c (below which there is no active phase) decreases with 

the range r; in fact, we find numerically that random process results into the action for a time slice. 


1 0.1612 

n c « - +- 

2 r 


(El) 


On figure [14] (solid lines) we show how for r = 6 the phase 
diagram is already close to the MF model (thin lines). Intu¬ 
itively, this is not surprising: When increasing the number of 
sites onto which an active site can depose its grains during a 
toppling event, these “neighbouring” sites are better and better 
described by a MF approach. 


It is also intuitively clear, that the same phenomenon takes 
place in high dimension d. Indeed, for d = oo, each active 
grain has 2 d neighbours, and supposing the state can locally 
be described by the number of grains and the activity, the as¬ 
sumptions in the derivation of MF theory become valid in the 
large-rf limit. 


Appendix F: Changing variables in the CSPI 


e -<5S[p*,p] _ / e -p*[p(t+5t)-p(t)] /gp*<5p 


= / e 


-p“[p(t+5t)-p(t)] 


1 + 


p{p)6t(e p * Sp -l) 


-1 +- 


(FI) 


Taking St small, and concatenating many such slices yields 
the path integral j T>[p]V[p*]e~ s ^ p - p \ with action 


S[p*,p] = Jp*(t)dtp(t) - H p [p*(t),p(t)] , (F2) 

U P {p*,p\=p(p)(e p * sp ~ l) . (F3) 

Consider now the examples discussed in section B There 
was pair annihilation with probability p(p) = ^p{p — 1) and 
Sp = —2; decay with rate pp and Sp = —1; finally particle 
creation (branching) with rate up and Sp = 1. Promoting the 
fields to operators, p —> p, and p* —> jy, the corresponding 
Hamiltonian reads 


In Ref. m the authors propose an alternative exact path- 
integral representation for stochastic systems, and show its 
equivalence to the CSPI formalism 6 . Their idea is to start 
from a rate equation for the particle number p on a site. 
Suppose in a time step St the particle number changes with 
probability p(p)St by a fixed amount Sp. To write down the 
path integral for a time slice St, the change p(t + St) = 
p(t) + Sp is enforced by an integral over the auxiliary field 
p*, f p , e~ p *^ t+st ^~ p ^~ sp y Averaging this quantity over the 


U p = V -\& 2pt -l]p(/5- 1) + p[e p+ —l]p + «[e pt —l]p . 

(F4) 

We remind that as in section urn the Hamiltonian is written 
in normal-ordered form, which is required in order to pass 
from the path integral to the operator formalism, and back. 
By construction, the commutation relations are canonical, 

[p, P f ] = I- (F5) 


6 The reader of Ref. QT) may benefit from the following remarks: (i) The 
sum in Eq. (2) contains both the pairs ( i , j) and (j, i). The unit vector ap¬ 
points from j to i. (iii) there is no factor of 2 in the definition of 7 after 
Eq. (3). (iv) In Eq. ( 6 ) ^ and v are exchanged. 


Let us compare ( |F4[ > to our Hamiltonian (151 for the same 
processes, reproduced here for convenience. 






-(fit) 


tf2-2 


d 


(fit) 2 fi — fitfi . 

(F6) 
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Demanding that Pi = H p uniquely identifies, already for v = 
0 , the operator identities 


a ^ 


a = e 


P ■ 


(FI) 


A few checks are in order: First, the commutation relations 
[a, a)} = 1 and [p,pF] = 1 are compatible. Since the trans¬ 
formation ( |F7[ ) is canonical, i.e. has Jacobian one, the path 
integrals also map onto each other. 

Finally, the term proportional to // is correctly reproduced; 
it is the only term which requires the commutation relations, 
and for which normal-ordering is important. 


We finally note that the Hamiltonian (46 1 for hopping, in 
the continuous limit, % = -/JVa'Va, maps onto 


U p = J -D\/^\/p + D{\/p^) 2 p. 


(F8) 


This expression can directly be derived from the equivalent 
of Eq. ( |F1| > for diffusion fl7) . Decoupling the term quadratic 
in p' with a noise yields Eq. (157 i. That these two equations 
are equivalent is non-trivial: the formalism used to derive 
Eq. ( |157| > works with coarse-grained fields, whereas p above 
is a microscopic field. The reason for this equivalence is that 
when coarse-graining Eq. ( |157[ ), its linear structure allows to 
pass coarse-graining through to the fields. 
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